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SUMMARY 


The objective of this research is to reduce the end-point vibration of a large, 
teleoperated manipulator while preserving the usefulness of the system motion. A master 
arm is designed to measure desired joint angles as the user specifies a desired tip motion. 
The desired joint angles from the master arm are the inputs to an adaptive P.D. control 
algorithm that positions the end-point of the manipulator. As the user moves the tip of 
the master, the robot will vibrate at its natural frequencies which makes it difficult to 
position the end-point. To eliminate the tip vibration during teleoperated motions, an 
input shaping method developed by Singer and Seering from MIT is presented. 

The input shaping method transforms each sample of the desired input into a new 
set of impulses that do not excite the system resonances. The method is explained using 
the equation of motion for a simple, second-order system. The impulse response of such 
a system is derived and the constraint equations for vibrationless motion are presented. 
To evaluate the robustness of the method, a different residual vibration equation from 
Singer’s is derived that more accurately represents the input shaping technique. The 
input shaping method will be shown to actually increase the residual vibration in certain 
situations when the system parameters are not accurately specified. Finally, the 
implementation of the input shaping method to a system with varying parameters is 
shown to induce a vibration into the system. To eliminate this vibration, a modified 
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command shaping technique is developed. The ability of the modified command shaping 
method to reduce vibration at the system resonances is tested by varying input 
perturbations to trajectories in a range of possible user inputs. By comparing the 
frequency responses of the transverse acceleration at the end-point of the manipulator, 
the modified method is compared to the original P.D. routine. The control scheme that 
produces the smaller magnitude of resonant vibration at the first natural frequency is 
considered the more effective control method. 


xiv 


CHAPTER I 


INTRODUCTION 

1.1 Motivation for this Research 

The industrial applications for robots with a relatively large workspace has 
increased significantly over the past few years. Most of the attention has been focused 
on the assembly of large space structures, the welding of airplanes and automobiles, the 
weaving of new composite structures such as submarines and, more recently, the 
inspection and removal of hazardous waste [9,22,23,48]. With an expanded workspace, 
the robot is often required to move large distances in a relatively short amount of time. 
The demand for a high speed robot with a large workspace usually requires long, 
lightweight links which are inherently flexible. This flexibility allows the link to store 
potential energy which is often returned to the system in the form of kinetic energy. 
Therefore, the end-point vibration of the manipulator, as well as uncertainty in the end- 
point position, can be directly related to the inherent flexibility of the links. 

For space and hazar dous waste environments, a remote operator is needed to 
perform most of the required tasks. The user often dictates the motion of the robot 
through the workspace and then precisely positions the end-point of the manipulator to 
perform a task. A teleoperated system is ideal for the user interaction requirements. 
Through the use of a joystick or other input device, the operator can specify the desired 
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trajectory in either end-point coordinates or joint coordinates. Therefore, the user’s 
desired motion can be transformed into actual robot motion. 

This research addresses the end-point vibration of a large, teleoperated 
manipulator. Present methods for reducing end-point vibration are discussed to 
determine the appropriate strategy to reduce tip vibration of the manipulator. 

1.2 Previous Methods for Reducing Vibration 

The problem of reducing end-point vibration is not a new one. Many different 
methods, both passive and active, have been investigated to eliminate unwanted 
oscillations. The most crude passive approach to eliminate vibration is to simply wait 
for the vibrations to stop after a desired motion. NASA originally used this method on 
their Space Shuttle Remote Manipulator System but found it to be costly in completion 
time requirements. Alberts and Book [1,2] experimented with a thin film of visco-elastic 
material that is applied to the structure surface and then covered with a very stiff 
constraining layer. When the beam is deformed, the visco-elastic material is sheared and 
thus energy is dissipated. The major drawback is that the vibrations are not eliminated 
but just reduced in amplitude. 

The majority of the strategies to eliminate end-point vibration involve active 
control structures. Different states of the system are measured and control efforts are 
based upon them. As with many effective control schemes, an accurate model of the 
system is required. Book [6] used distributed and lumped parameter models for each 
arm segment. The models are combined using homogeneous transformations and then 


numerical techniques are used to derive frequency domain information. Book [7] later 
derived the recursive dynamic equations for a flexible manipulator. However, the 
recursive method must be evaluated symbolically to obtain the dynamic equations in 
closed form. Finally, Book [8] reviews the mathematical representations commonly used 
in modeling flexible systems. 

Hastings and Book [17] extended active control methods by including strain 
feedback in the control structure. Their experiments showed that strain feedback can 
reduce the residual vibration during settling time. However, they concluded that the 
vibrations are inevitable with a feedback control scheme because the feedback control 
signal contains high frequency components, which excite the system resonances. 

Montgomery, Ghosh and Kenny [33] propose torque-wheel actuators to reduce 
overshoot in the Space Shuttle Remote Manipulator. Their method uses an inertial device 
to assist in reducing end-point vibration when following telerobotic commands. The 
results from their experiments indicate that the torque-wheel can produce a vibration of 
significant amplitude to diminish the original vibration while under teleoperated control. 
However, this procedure requires external devices to be mounted on the robot and is 
shown to work only for an abrupt stop command. 

Tewani, Walcott and Rouch [45] suggest using a dynamic absorber as a viable 
means for suppressing vibrations of a system. The method involves the combination of 
passive elements, active elements and an absorber mass to apply a controlling force to 
the system. Using a disturbance rejection control strategy, the amplitude of vibration 
was significantly reduced. However, this method would be difficult to implement on a 
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tcleoperated system because the disturbances are never known exactly. By implementing 
a Linear Quadratic Regulator (LQR) controller, a reduction in amplitude of vibration was 
still observed. 

Presently, Lee and Book [28,29] are studying the effects of inertial forces to 
suppress the vibration of a large, flexible robot. By mounting a small robot at the tip 
of the large robot, damping forces are generated to accommodate the inertial forces 
generated when the tip vibrates. Simulation results have shown the effectiveness of the 
damping forces and the controller designed. Using deflection rate control, vibrations are 
damped in half the time required with passive control. Currently, Lee is conducting 
experiments using the prescribed robot configuration to verify simulation results. 

Singer and Seering [39,40] presented a method of generating shaped command 
inputs to reduce end-point vibration. Unlike the previous methods which measure system 
states to reduce vibration, Singer’s method utilizes system information to alter input 
commands to the actuators. Each commanded impulse is appropriately distributed into 
a multiple impulse input whose characteristics are based on the system’s natural 
frequencies and damping ratio. This procedure, in effect, filters out frequency 
components near the system’s resonances to avoid exciting the system. Later, Singer and 
Seering [38,42] show that the input shaping idea is effective with teleoperator inputs. 
However, their system was limited to a beam operating in the horizontal plane so that 
the natural frequency was constant for a given experiment. 
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1.3 Characteristics of a Flexible System 

The equations of motion of a flexible system can be quite complex when 
compared to the rigid body counterpart. Lagrange’s equations of motion produce an 
infinite number of vibrational modes and frequencies for just a simple Euler beam 
[14,31,32]. The ability to regulate these modes is limited by the bandwidth of the 
control system. Therefore, only a finite number of modes of a flexible system can be 
controlled. 

Nonlinear effects due to large tip velocities may also need consideration. The 
centrifugal and Coriolis accelerations may generate additional tip vibration that is not 
controllable with conventional feedback schemes. The deflection of the manipulator due 
to vibration may also exhibit nonlinear properties if the amplitude of vibration is large 
relative to the link length. 

The nonminimum phase characteristic of flexible systems must also be considered 
for tracking teleoperated inputs. Kwon and Book [24,25] show that the transfer function 
between the input at a joint and the end-point position has zeros with positive real parts. 
These right-half complex plane zeros cannot be canceled using conventional feedback 
control algorithms. This prevents the feedback controller from having asymptotic 
tracking stability. 

Finally, the issue of noncolocated control must be mentioned. Noncolocated 
control occurs when the control effort and the sensing of the system states do not occur 
at the same point in the system. An example is when a robot is actuated at the joints but 
the end-point position is sensed at the tip. Often, nonminimum phase characteristics arise 


5 



when a system is noncolocated. Therefore, a noncolocated, state feedback control 
scheme must be designed with some consideration for system instability. 

1.4 Method of Approach 

The main focus of this research is the alleviation of end-point vibration in a 
telerobotic system. To avoid the complex derivation of dynamic equations, an adaptive 
control strategy derived by Yuan is used [49,50]. His controller was based on the error 
between the desired joint angles and the actual joint angles. By devising an operator 
input device that commands desired joint angles, Yuan’s controller was easily 
transformed to work, as a teleoperated control system. 

This control structure is then compared to a control scheme that performs input 
shaping to the joint error signal. First, the original input shaping method developed by 
Singer is implemented and shown to produce a vibration in the system. To prevent this 
vibration, a modified command shaping technique is developed. 

The advantages of using this new modified method over the original adaptive 
control scheme are then examined. By giving the two control methods identical input 
trajectories, the frequency responses of the transverse acceleration at the end-point of the 
manipulator are compared. The control scheme that produces the smaller magnitude of 
resonant vibration at the first natural frequency of the system is considered the more 
effective control method. 
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CHAPTER n 


EXPERIMENTAL TEST BED 


2.1 Description of the Robotic Manipulator 

The flexible manipulator used in the Flexible Automation Laboratory at Georgia 
Tech was designed by a Master’s student in 1986. Wilson [47] designed the robot to 
carry a payload of around 100 pounds, to reach second story windows and to be 
lightweight. To meet the design specifications, a two degree-of-freedom manipulator was 
built using 10 feet long, aluminum links. By choosing aluminum, the structure remains 
lightweight with good strength properties. 

The manipulator, named RALF (Robotic Arm, Large and Flexible), is shown in 
Figure 2.1. Both links are maneuvered using single ended hydraulic actuators. Actuator 
one is attached directly to a collar on link one while actuator two positions link two using 
a four-bar parallel mechanism. This parallel mechanism gives much needed support to 
allow the robot to lift the 100 pound payload requirement. The mechanism also stiffens 
the overall structure which raises the system natural frequencies. However, this parallel 
linkage produces nonlinear effects (e.g. dead bands) and couples the dynamics of the two 
links together. 

The nonlinear dynamic equations of motion were derived by Lee [26] using 
Lagrange’s equation. These nonlinear equations were evaluated in symbolic form and 
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then, using single value decomposition, the constrained dynamic system was solved. 



Figure 2. 1 Two-Link Flexible Manipulator 


Since hydraulic actuators are used, the joint motions are limited. Joint one has 
a range from 35° to 110°, which is measured from the horizontal plane. Joint two is 
limited to angles of 55° to 108° when measured relative to link one. Even with limited 
joint motion, the workspace, shown in Figure 2.2, is still quite large. Notice how the 
first link of RALF is able to pass through the vertical axis, y ot that intersects joint one. 
This allows the center of gravity of the manipulator to move from one side of the axis 
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to another which can create unusual dynamics for the actuator connected to link one. 
Since the forces required for a given motion are different in each region, the actuator 
dynamics play an important role in the frequency analysis conducted in Chapter IV. 



Now that the physical structure of RALF has been discussed, the forward and 
inverse kinematics are presented. The derivation of the kinematic equations is best 
handled using matrix transformations. By representing each degree of freedom of the 
manipulator with an individual coordinate frame, a matrix transformation is created 
[1 1 ,44], Each matrix transformation relates the current coordinate system to the previous 
one. The overall transformation, relating the end-point to the base of the robot, is 
obtained by multiplying the individual matrix transformations together. For a flexible 
system, Book [7] demonstrated that the overall matrix transformation is actually the 
superposition of a rigid body transformation and a flexible transformation. TTus stems 
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from the fact that the position of a point in a flexible system can be described using rigid 
body coordinates along with the modes of vibration. The rigid body coordinate frames 
used for RALF are shown in Figure 2.3. 



For slow motions, the flexible modes of vibration may not be excited. In this 
situation, the rigid body transformations give a reasonable approximation of the tip 
position. Since modeling of the flexible dynamics is not the main emphasis of this 
research, only the rigid body transformations are considered. The rigid body 
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transformations for RALF are given in Appendix A. 

The forward kinematics problem can be stated as follows: given the joint variables 
of the robot, determine the position and orientation of the end-effector. For RALF, only 
the (x,y) coordinates of the end-point are important. From the overall matrix 
transformation given in Appendix A, the last column yields the desired kinematic 
equations. The forward kinematic equations are 


x * ljcos(0, +0 2 ) + L,cos(0,) - I3 sin(0j) (2.1) 

y « L z sm(Q l + 0 2 ) + LjSinfOj) + LjCosfO,) . (2.2) 

The inverse kinematics problem is more difficult than the forward kinematics 
solution. The inverse situation simply stated is: given a desired position and orientation 
for the end-effector, determine all the possible joint configurations that achieve the 
desired position and orientation. Since the resulting kinematic equations are nonlinear, 
there is no guarantee of finding a unique solution or even finding a real solution at all. 

The existence of a solution defines the workspace of the manipulator. The lack 
of a solution indicates that the desired position and orientation are not within the robot’s 
workspace. 
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From Appendix A, the inverse kinematic equations are 




atan2 


kl x -i,/ 4^0? 


(2.3) 


0j * atan2 


' -x(L I sin(6 2 ) +£3) ♦yfLgCosfS^) *L t ) \ 
k xCI^co^e^+I^+yCIjSmCej)*!,) J 


(2.4) 


where 


K 


— 


- 1} ~ Is- 


(2.5) 


2.2 Design of an Input Device 

The main function of the input device in teleoperation is to assist the operator in 
accurately maneuvering the end-point of the manipulator. By sensing the user’s desired 
motions, the device should convert these motion commands to movements of the end- 
point. The design of such an input device should be functional and suitable for the 
specific application. 

Fischer, Daniel and Siva [13] discuss many guidelines for the design of input 
devices for use in teleoperation. The engineering group emphasized the necessity of the 
feedback of key information such as position and forces to avoid damage to the robot. 
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Since bracing of the manipulator is not considered and the workspace of RALF is 
unobstructed, only the relative position of the tip in the workspace is necessary. 

Since RALF has only two degrees of freedom, the design of the input device was 
quite simple. A two link scaled model of RALF was designed with linear potentiometers 
placed at each joint to record the user joint commands. By moving the master arm, the 
operator can think of positioning the end-point of the manipulator in Cartesian 
coordinates while the input device functions in joint coordinates. This strategy of 
measuring joint commands directly works well with the adaptive control routine 
developed by Yuan discussed in Chapter I. 

The overall teleoperated system designed to position the end-point of RALF is 
shown in Figure 2.4. The workspace of RALF is mounted as a backdrop to give the 
operator a scaled picture of the allowable workspace. This arrangement makes it easier 
for the user to make relative maneuvers within the workspace. 

2.3 Interfacing the Telerobotic System 

The coordination of the master arm with the slave arm in a telerobotic system is 
important to produce desired results. Fiala [12] describes a logical architecture for 
connecting teleoperation input devices to the telerobotic control hierarchy. A method of 
handling control information is presented which allows many system components to 
access the information simultaneously. He also explains the two main classes of 
teleoperation input devices, joint-space devices and Cartesian devices. 
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The effects of varying system parameters on the ability of the operator to position 
the end-point of the manipulator should be understood. Hannema and Book [15] discuss 
moving the end-point of a manipulator from one point to within a certain tolerance band 
surrounding a desired point. They examined the effects of backlash, Coulomb friction 
and bandwidth on the ability of an operator to position the end-point of a manipulator. 
They showed that a linear model, relating the task parameters of distance and width to 
performance in task completion time, could be made. Their experimental results 
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revealed an improvement in task completion time as the joint bandwidths were increased 
from 1 to 3.5 Hz. 

Uebel, Ali and Minis [46] also investigated the effects of bandwidth on operator 
performance using a Robotics Research Corporation slave arm with a Kraft Telerobotics 
master arm. Varying the joint bandwidths from 0.5 to 2 Hz, experienced operators 
performed five repetitions of a peg-in-a-hole task for three different bandwidths. Their 
results also show a decline in operator performance as the bandwidth is decreased. 

The final issue that must be considered when interfacing telerobots is the effects 
of time-delay. Niemeyer and Slotine [34] address the problem of time-delays in 
telerobotics by using an adaptive control strategy. By using an adaptive controller, the 
manipulator bandwidth is not limited to the bandwidth imposed by the transmission 
delays. However, the type of data that is transmitted (e.g. forces, torques) can alter the 
behavior and stability of the overall system. 

2.4 System Hardware 

The telerobotic system used for experiments at Georgia Tech consists of the input 
device discussed in Section 2.2 functioning as the master, and RALF as the slave. The 
main control unit is a MicroVAX II made by Digital Equipment Corporation. The 
necessary connections to the sensors, hydraulic valves, amplifiers and master arm are 
made using an analog/digital (A/D) board made by Data Translation (see Appendix C for 
equipment list). The A/D boards can sample a single channel at 6000 Hz. However, 
eight channels sample at only 300 Hz. A second A/D board is available but a one 
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milliseco nd delay is required to switch from one board to the next board. 

After the computation time for the control routines is considered, the sampling 
frequency is reduced to 50 Hz. The frequencies of RALF to be controlled range from 
3.7 to 5.5 Hz. From Nyquist criterion, the A/D sampling rate is still sufficient to 
control these system frequencies. 

The adaptive control algorithms, the modified command shaping routines and the 
A/D routines were written in VAX C. Previous control routines were written in 
FORTRAN (FORmula TRANslation) and had sampling rates of around 150 Hz. 
Obviously, the FORTRAN code runs almost three times faster than the C code. Test 
programs verify that the Micro VAX n machine is optimized to run Fortran code. The 
VAX C language was chosen for portability to personal computers (PC) where future 
control will be implemented. 


CHAPTER m 


MODIFIED COMMAND SHAPING 


3.1 Historical Perspective 

The original command shaping technique to reduce system vibration was 
introduced by Singer and Seering in 1988 [39]. The procedure involved modifying 
desired system inputs so that a system completes the requested motion without vibration. 
The method was tested on a computer model of the Space Shuttle Remote Manipulator 
System (RMS) developed by Draper Laboratories. Simulation results from a variety of 
commanded moves suggested that significant vibration reduction can be obtained using 
the method. Subsequent documents [40,41] verified the ability of the method to reduce 
end-point vibrations. 

By 1989, Singer showed that the procedure worked with multiple mode systems 
and was effective on telerobotic systems as well [38,42]. Singhose, Seering and Singer 
expanded the method using a vector diagram approach to determine the appropriate 
shaping strategy given some allowable residual vibration amplitude [43]. Hyde and 
Seering [21] extended these results to the solution of a group of simultaneous non-linear 
impulse constraint equations. 

Recently, the input shaping method was implemented on an overhead crane to 
reduce oscillatory motion of the object being moved. Noakes and Jansen [35] 
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gene ralized the theory associated with oscillation-damped trajectories to a system with 
simply suspended loads. Previously, a crane operator moved the suspended object slowly 
and allocated time for pendulum oscillations to damp out. After implementing the 
oscillation-damping algorithms, they were able to position 55-gallon drums in a U-shaped 
path with insignificant oscillations. 

3.2 Derivation of Constraint Equations 

The original input shaping method involves the manipulation of a desired input 
command. Each sample of the input command is replaced by a sequence of impulses that 
do not excite the system natural frequencies. Knowing the impulse response of the 
system, constraint equations can be derived that yield the appropriate amplitudes and 
starting times of the impulse sequence. 

The constraint equations can be derived from the impulse response of a simple 
linear, time-invariant second-order system. Consider the spring-mass-damper system 



Figure 3.1 Second-Order System 
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shown in Figure 3.1. The vibratory response of this system to an impulse input is 


. -{«.(<-»#) 

x(f) * — r — smfw.Vl "C 2 (f -*<>)) 
t/W 5 


(3.1) 


where A is the amplitude of the impulse, w. is the natural frequency of the system, f is 
the damping ratio of the system, t is the time and to is the time when the impulse occurs. 
This result is derived in Appendix D. Using Equation (3.1), the position response, x(r), 
for the second-order system is specified for time, t ^ t 0 . 

If the system is given a two impulse input, the vibratory response is 

x(r) * BjSinfaf +<J> A ) +fl 2 sin(af +^ 2 ) (3-2) 


where 

» . <*-«■•> (3.3) 

/W 5 


a * Uj 7 1 - C 2 


( 3 . 4 ) 


* - “..'otv'l -C 2 . 


( 3 . 5 ) 
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The two impulse response given in Equation (3.2) can be simplified to yield 




(3.6) 


where 


B ^ +B 2 cos(4> 2 )] 2 +[B,sin(<|> 1 ) +5 2 sm(4> 2 )] 2 


(3.7) 


^ = tan -1 


k BjCos(4>j) + B 2 cos(<}> 2 ) J 


(3.8) 


The summation of two sinusoids is proven in Appendix E. 

Since the system is linear and time-invariant, the results from Equations (3.7) and 
(3.8) can be generalized to the N impulse input case. The amplitude and phase of 
vibration for the N impulse input case are 


B 


amp 


N 


J B k cos(<t> t ) 


12 U t sin(<f> t ) 


(3.9) 


/ N 

52 £ 4 sm(4> 4 ) 




-l 




^ 52 fl 4 cos(<t> 4 ) j 


(3.10) 
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Since the purpose of the input shaping method is to eliminate vibration, the 
amplitude of vibration, Equation (3.9), must equal zero after the last impulse occurs. 
This only happens if both the squared terms are independently zero since the sine and 
cosine functions are linearly independent. The resulting equations are 

B x cos(4» x ) + £ 2 cos(4> 2 ) +... + B n cb&($ n ) - 0 (3.11) 

£ 1 sin(<t> 1 ) + B 2 sm(4> 2 ) ♦••• + ® (3-12) 


including 






(3.13) 


and 




(3.14) 


where A k is the amplitude of the it* impulse, r * is the time when the final impulse 
occurs and t a is the time when the k* impulse occurs. 


21 



Substituting the equations for B k and 4> k into Equations (3.11) and (3.12), the 


constraint equations become 


£ 008(0^1 -C 2 t 0k ) « 0 (315) 


f;^e' C “* ( ‘"' W sin((o |i v'l -C 2 / 0 *) - (316) 

*«i 


For the constraint equations to produce the correct impulse sequence to eliminate 
vibration, the natural frequency and damping ratio of the system must be known exactly. 
Since these system characteristics are not precisely known, their robustness is included 
as a constraint. The robustness constraint with respect to natural frequency is found by 
taking the partial derivative of Equations (3.15) and (3.16) with respect to and setting 
the result equal to zero. Likewise, the robustness constraint with respect to damping 
ratio is found by taking the partial derivative with respect to f and setting this result 
equal to zero. After performing the described differentiation, the resulting constraint 
equations are the same. Therefore, setting the partial derivative with respect to w, equal 
to zero is equivalent to setting the partial derivative with respect to f equal to zero [38]. 
These new constraint equations form the first-derivative robustness criterion. 
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Mathematically, the equations for the first-derivative robustness criterion are 


E Vo* ‘ W «* ( «„,/ 1 - C 1 »«) -0 (3.17) 

k-l 


« 0 . < 318) 


Higher derivative constraints are obtained by differentiating Equations (3.15) and 
(3.16) to the desired order. The m ^-derivative robustness constraint equations are 


E 4k ( f o * )* e C cos( w.V 1 - c 2 t 0k ) « 0 

k-l 


(3.19) 


E ^ ( <o« ‘ ‘ w “>( “./w 3 »„) - o. 

t-l 


(3.20) 


The length of the impulse sequence is now determined by the number of 
unknowns in a given set of constraint equations. For any given set, there will always be 
two more unknowns than equations. To alleviate this dilemma, the starting time of the 
first impulse is arbitrarily chosen to be time zero and the amplitudes are normalized so 
that they sum to one. This particular normalization ensures that the overall amplitude 
of the new impulse sequence is the same as the amplitude of the desired input command. 
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3.3 Calculation of Impulse Amplitudes 

Now that the robustness constraint equations have been determined, the impulse 
amplitudes and starting times can be solved. For the two impulse input, the zero *- 
derivative constraint equations are utilized, which are 


BjCosC^) + B 2 cos(4> 2 ) ■ 0 (3.21) 


^sinW,) + B 2 sin(<|> 2 ) - 0 


(3.22) 


and 


B , Ak(i3 * 


/W 1 


(3.23) 


where A k is the amplitude of the k * impulse, t N is the time when the final impulse 
occurs (i.e. and t 01 is the time when the first impulse occurs. 

Since any equation involving sines and cosines is transcendental, there are an 
infinite number of possible solutions to Equations (3.21) and (3.22). Therefore, only the 
solution that yields the shortest time duration and a positive amplitude for all the 
impulses is chosen. For the two impulse input case, a complete derivation of the solution 
process is furnished in Appendix F. 
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The resulting solution for the two impulse case is 


A. - — (3.24) t 01 - 0 (3.25) 

1 + Af 


** - 


M 

1 +M 


(3.26) 


*02 


n 

o./w 5 


(3.27) 


where 


M • e 



(3.28) 


The ability of the input shaping method to eliminate vibration can be demonstrated 
graphically. Consider the input in Figure 3.2 whose characteristics are given in 
Equations (3.24) through (3.27). 



Figure 3.2 Two Impulse Input 
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The second-order system response to each of the individual impulses in Figure 3.2 is 
shown in Figure 3.3. 



Figure 3.3 System Response to Each Impulse 


Since the system is linear and time-invariant, a linear combination of two inputs 
results in a response that is a linear combination of the two responses. Therefore, the 
net system response to the two impulse input is shown in Figure 3.4. Since the natural 
frequency and damping ratio of the system are exact, there is no vibration of the system 
after the second impulse. The effects of parameter uncertainty on the amount of residual 
vibration is discussed in Section 3.4. 
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0.8 


Overall System Response to Two Impulse Input 

w • « a 




Figure 3.5 Single Impulse Response vs. Two Impulse Response 


27 





Figure 3.5 shows a comparison of the impulse response of a second-order system 
to the two impulse input response. Since the inputs have the same amplitude, i.e. they 
sum to one, the responses can be compared to determine which is more desirable. For 
this simple case, the two impulse input is preferred since it eliminates the vibration after 
the second impulse occurs. The ability to completely eliminate the vibration is attributed 
to an ideal second-order system model which provides exact system parameters. 

To solve the three impulse input case, the zero* and first derivative constraint 
equations are evaluated. The complete derivation for the three impulse case is in 
Appendix G. The solution for this case is 




1 

1 + 2Af ♦ M 2 


(3.29) 


t 0l = 0 (3.30) 


A = — 

2 1 + 2Af + M 2 


(3.31) 


M 2 


A, = - (3.33) 

1 + 2M +M 2 


*02 = 


1 - c 2 


(3.32) 


2 ic 


*03 


i - c 2 


(3.34) 
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Figure 3.6 shows the three impulse input case whose characteristics are given in 
Equations (3.29) through (3.34). 



Finally, the four impulse input case is solved. To obtain the amplitudes and 
starting times of the four impulses, the zero *, first and second derivative constraint 
equations must be solved. Appendix H contains the complete derivation for this case. 
The results from this appendix are 


A, = !— J (3.35) 

1 -t3Af ♦3M 2 »M 5 


• 0 (3.36) 


A . - (3.37) 

1 +3Jf + 3M 2 +Af 3 


* 02 


It 

o./w 5 


(3.38) 
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3M 2 

1 + 3M + 3Af 2 +M 3 


2 * 


(3.40) 


^3 


(3.39) 


w./W 2 


; (3-41) 

1 + 3Af + 3A/ 2 + Af 5 


*04 


3n 

«.Vi -c 2 


(3.42) 


Figure 3.7 shows the four impulse input whose characteristics are given in Equations 
(3.35) through (3.42). 
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3.4 Robustness of Constraint Equations 

In Section 3.3, the two impulse input completely eliminated the vibration when 
the natural frequency and damping ratio of the system were known exactly. For most 
physical systems, the exact parameters are seldom known. Thus, there is some residual 
vibration after the last impulse has occurred. 

To determine the amount of residual vibration, a vibration error expression must 
be defined. The error, denoted err, is written as the ratio of the actual multiple impulse 
response magnitude to the actual impulse response magnitude of the second-order system. 
The error expression is defined only for time after the multiple impulse input has 
occurred to ensure that the system has received identical amplitude inputs. 
Mathematically, the vibration error is written as 


err 


|*.CO| 


for t a t ok 


(3.43) 


where k is the number of impulses. The residual vibration is just the vibration error 
expressed as a percentage. 

The deviation in the actual system parameters from the design parameters can now 
be quantified using Equation (3.43). The actual system response to a multiple impulse 
input can be computed and related to the actual impulse response of the second-order 
system. By studying the deviations in the natural frequency and damping ratio from the 
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design parameters, their effects on the vibration error can be better understood. 

This new definition of the vibration error expression is much different than the 
one originally stated by Singer. He defined the vibration error expression as "the 
maximum amplitude of the residual vibration during a move as a percentage of the 
amplitude of the rigid body motion." This definition is expressed mathematically with 
Equation (3.9) divided by the sum of all the A k [39], which is always unity. The main 
problem with Singer’s definition is that it does not accurately represent the ability of the 
input shaping method to reduce vibration. In some instances, the input shaping method 
can actually increase the residual vibration of a system. 

To prove this point, the deviation in actual natural frequency, u a , from the design 
natural frequency, «„ is analyzed. The vibration error defined by Singer for the two 
impulse input case is shown in Figure 3.8. Singer states that an acceptable vibration 
error level is less that 5% residual vibration for a second-order system. Therefore, the 
two impulse input is robust for a frequency variation of less than 5 % . From Figure 3.8, 
the residual vibration curves decrease in magnitude for an increase in the value of 
damping ratio, f. The graph also shows that the magnitude never exceeds 100% and 
therefore the input shaping method can never increase the residual vibration of a system. 
This implies that the input shaping method reduces residual vibration for any variation 
in natural frequency. However, the new definition of residual vibration clearly shows 


this is not the case. 


The vibration error from Equation (3.43) for the two impulse input case is 



1 


|*«(0| 

Af + 1 \ 

K J 


(3.44) 


where is the actual natural frequency of the system and u„ is the design natural 
frequency of the system. Equation (3.44) is derived in Appendix I. Figure 3.9 shows 
the new vibration error as a function of normalized frequency, caju,,. 


Vibration Error vs. Normalized Frequency 
N- ? 



Figure 3.8 Vibration Error vs. Normalized Frequency 
Singer’s Definition for Two Impulse Input 
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Vibrotion Error vs. Normalized Frequency 
N - 2 



Figure 3.9 Vibration Error vs. Normalized Frequency 
Two Impulse Input 


The residual vibration now increases for an increase in damping ratio, f. This 
fact may seem incorrect since the overshoot of a second-order system increases with a 
decrease in damping ratio for a step input. However, analysis of the impulse response, 
given in Equation (3.1), verifies the results displayed in Figure 3.9. 

Figure 3.9 also shows that replacing the original impulse by a sequence of 
impulses can actually have a negative effect if the actual natural frequency is over 1.5 
times the design frequency. Figure 3.10 shows a second-order system response to a 
single impulse input compared to a two impulse input when the actual natural frequency 
is twice the design frequency. The residual vibration after the second impulse is actually 
worse than if the system had only been given the single impulse. Therefore, input 
shaping can have a detrimental effect when large errors in natural frequency are present. 
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Single Impulse Response vs. Two Impulse Response 

tt, - 2*«, 



Figure 3. 10 Single Impulse Response vs. Two Impulse Response 


Since the two impulse input case is only robust for deviations in natural frequency 
of less that ± 5%, the robustness of the three and four impulse cases is of interest. 
Figure 3.11 shows the vibration error versus normalized frequency for the three impulse 
input. For this case, the input is robust for deviations in natural frequency near ±10%. 
However, the ability of the input shaping method to produce detrimental effects is much 
more pronounced. 

By separating the input into four new impulses, Figure 3.12 shows that the 
method is robust for deviations in natural frequency close to ± 20%. Nevertheless, the 
possible adverse effects are even more noticeable. If the actual frequency is more than 
40% larger than the design frequency (for f = 0.2), then the input shaping method is of 
no practical use. 


35 





Figure 3.11 Vibration Error vs. Normalized Frequency 
Three Impulse Input 



Figure 3.12 Vibration Error vs. Normalized Frequency 
Four Impulse Input 
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The input shaping method must also be robust for deviations in damping ratio. 
Section 3.2 stated that the derivative of the constraint equations with respect to damping 
ratio yields the same constraint equations as the derivative with respect to natural 
frequency. Therefore, robustness in damping ratio is already accomplished when 
robustness in natural frequency has been considered. Figure 3.13 shows the vibration 
error versus normalized damping ratio for the multiple impulse cases. 



Vibrotion Error vs. Normolized Domping Ratio 
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Figure 3.13 Vibration Error vs. Normalized Damping Ratio 

To evaluate the error expressions for the deviation in damping ratio of a second- 
order system (see Appendix J), a specific value of damping ratio is required. For this 
example, the design damping ratio is 0.05. For all three input cases, large deviations in 
damping ratio do not have a significant effect on the residual vibration. This fact is 
comforting since the damping ratio for a complex system may be hard to measure. 
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3.5 Position Dependent System Parameters 

The robustness of the constraint equations, discussed in Section 3.4, demonstrated 
the ability of the input shaping method to reduce vibrations even with deviations from 
the design system parameters. However, Singer’s original input shaping technique does 
not address the issue of changing system parameters. For RALF, the natural frequency 
and damping ratio are functions of position, i.e. joint angles. Therefore, a modified 
command shaping technique was developed to accommodate varying system parameters. 

First, the implementation of the input shaping technique to a discrete-time system 
is presented. Figure 3.14 shows a simple block diagram of the input shaping method. 
For each sample of the input, N output impulses are generated. From Section 3.3, the 



Figure 3.14 Input Shaping Block Diagram 


time period between output impulses is the same. This time period, denoted delT, is 


delT 




(3.45) 


which is a function of both the natural frequency, &>„, and the damping ratio, f. 
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To u tilize this time period information in a discrete-time system, the continuous- 
time data must be represented in discrete-time. From discrete-time signal processing, 
a continuous-time signal, x(r), is represented mathematically as a sequence of numbers, 
x[n], where n is strictly an integer. To transform the continuous-time period delT into 
a discrete-time period dein, the sampling rate of the discrete-time system, f„ is used. 
The equation to perform this transformation is 

debt - int ( delT */,) (3-46) 

where the int function truncates the argument to an integer. 

For the input shaping method, the discrete-time period, deln, never changes 
because the system parameters are assumed constant. But when the input shaping method 
is applied to a system that has time varying parameters, the continuous-time period, delT, 
becomes time varying as well. A significant change in delT will result in a change in the 
discrete-time period, deln, which produces an undesirable vibration in the system. The 
amount of change in the continuous period that causes this change in the discrete period 
is a function of sampling rate since deln is strictly an integer. 

For example, a four impulse output, shaping scheme is applied to a system that 
causes a change in deln from four to five. When each input sample is shaped into four 
output impulses, the method would produce a steady-state impulse output shown in 
Figure 3.15. Each output impulse is designated {a.b} where a indicates the discrete-time 
location of the input sample responsible for the four output impulses and b indexes the 
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four resulting output impulses. After examining Figure 3.15, it is obvious that the 
change in deln has caused gaps in the output for discrete values of n. At n= 4, for 
example, only three impulses are contributing to the overall output. To make matters 
worse, this problem is repeated five more times at a discrete-time period near the 
system’s natural period. This phenomenon induces a vibration into the system that is 
caused solely by the application of the input shaping method to a system with time 
varying parameters. 

This induced vibration is also present when the value of deln decreases. Consider 
a four impulse output, shaping scheme that is applied to a system that causes a change 
in deln from five to four. The resulting steady-state impulse output is shown in Figure 
3.16. For this situation, a surplus of output impulses is generated at a discrete-time 
period near the system’s natural period. These extra impulses also cause a vibration that 
is produced by the input shaping method. 

To eliminate the induced vibration, a modified command shaping method is 
proposed to make the impulse output more uniform when a change in deln in 
encountered. To compensate for a change in the discrete-time period, extra impulses are 
added for an increase in deln and impulses are removed for a decrease in deln. The 
choice of which impulses are affected is based on the number of output impulses from 
the shaping algorithm and the old and new values of the discrete-time period. 

The modified command shaping method can be explained by designating the 
discrete-time value when the discrete period increases as n=0. For the next N - 1 samples 
of the input, i.e. 0 < n £ N- 2, the modified command shaping technique shapes each 
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sam ple using both the old 2 nd new V2lues of deln to create 2 smooth steady-state impulse 
output. Using the new value of deln , the input sample is shaped to create N output 
impulses that are added to the overall output at their respective discrcte*time values. 
Using the old value of deln, the same input sample is also shaped to create N output 
impulses. However, only the last N-(n+ 1) output impulses are added to the steady-state 
output at their respective discrete-time values. For discrete-time values of n ^ N- 1, 
each sample of the input is shaped normally using the new value of deln to generate the 
N output impulses. 

The modified command shaping method also works for a decrease in the discrete- 
time period, deln. For this situation, the input sample is shaped only once using the new 
value of deln to produce the four output impulses. Instead of adding all four of the 
output impulses, only the first (n+ 1) output impulses are added to the steady-state output 
at their respective discrete-time values. By manipulating the overall output in this way , 
the extra impulses that are added for the case when deln increases are the same impulses 
that are removed when deln decreases. 

One final case to consider is when the value of deln changes more than once 
within one discrete-time period. For this situation, a new modified technique must be 
devised. For instance, if the discrete-time period length changes from one value to 
another and back again, the best method to smooth the steady-state output may be to 
ignore the change in discrete period if it is relatively short. 

To understand the modified command shaping procedure, consider the example 
given in Figure 3.15. The value of deln increases from four to five for this input 
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shaping scheme that produces four output impulses. Since N = 4 for this case, the next 
three (i.e. N-l) input samples will be shaped twice. At discrete-time n=0, the input 
sample is shaped using the new value of deln (i.e. 5) to create four (i.e. N ) output 
impulses that are added to the overall output. At the same discrete-time value, the input 
sample is shaped using the old value of deln (i.e. 4) to create four (i.e. N ) output 
impulses. However, only the last three (i.e. iV-(n+l)) impulses are added to the overall 
output at their respective discrete-time values. For the next discrete-time value, i.e. 
71 = 1 , the input sample is shaped using the new value of deln to create the usual four 
output impulses that are added to the overall output. When the input sample is shaped 
using the old value of deln, only the last two (i.e. N-(n+ 1)) output impulses are added 
to the general output. This process of shaping the input samples twice is repeated until 
the discrete-time value, n, is greater that N- 2. After n > N- 2, the shaping continues 
normally using only the new value of deln to produce the output impulses. 

This modified command shaping technique is demonstrated on the two examples 
discussed previously when the number of output impulses is four. Figure 3.17 shows the 
modified shaping technique implemented for the case when the discrete-time period 
increases from four to five. The impulses due to the modified command shaping 
technique are darkened to show emphasis only. Figure 3. 18 shows the modified shaping 
technique applied to the case when the period decreases from five to four. The impulses 
that are created but not added are drawn in the figure without tails to distinguish them 
from the normal impulses. The success of this modified input shaping technique is 
discussed in Chapter V. 
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Figure 3.15 Steady-State Output for an Increase in Discrete-Time 
Period Using Input Shaping 
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Figure 3.16 Steady-State Output for a Decrease in Discrete-Time 
Period Using Input Shaping 
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Figure 3.17 Steady-State Output for an Increase in Discrete-Time 
Period Using Modified Command Shaping 
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Figure 3.18 Steady-State Output for a Decrease in Discrete-Time 
Period Using Modified Command Shaping 
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CHAPTER IV 


DYNAMIC ANALYSIS OF RALF 


4.1 Previous Frequency Analysis 

The vibration analysis of RALF has been described by several investigators in the 
last five years. Huggins and Lee have performed the most research on verification of 
modeling methods and control algorithms for RALF. Huggins [19,20] conducted 
extensive finite element analysis to verify mathematical and experimental models used 
to determine the system modes and natural frequencies. Lee [26,27] derived a nonlinear 
model of RALF and then verified it through simulations using TREETOPS, a computer 
software package. However, these investigators only conducted research on a limited 
workspace of RALF. For a teleoperated system, the frequency data for the whole 
workspace is needed to implement the modified command shaping technique. 

4.2 Experimental Setup and Procedures 

The experimental determination of system resonances and damping ratios using 
digital Fourier analyzers is well published. Ramsey [36,37] discusses the importance of 
understanding the dynamic behavior of vibrating systems. Using Fourier analyzers, he 
explains many effective measurement techniques for modal analysis of different vibrating 
systems. The frequency response of a second-order system determined from several 
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different test inputs is also presented and the trade-offs for each method is discussed. 

Hewlett-Packard has published many application notes about using their Digital 
Signal Analyzer (HP3562A) for modal analysis [3,4,5]. After reviewing these 
documents, random noise was chosen as the input to determine the desired system 
properties. Using the P.D. feedback control algorithm running on a MicroVAX to hold 
RALF is a desired joint configuration, the random noise signal from the analyzer was 
added to the control signal to stimulate the system. The system was excited at the 
control level instead of using an external shaker to include actuator dynamics that might 
influence the results. The circuit that combined the two signals is shown in Appendix 
K for reference. An accelerometer, mounted at the tip of RALF, measured the 
transverse acceleration of the second link which was then recorded by the analyzer. 

For a given joint configuration, the power spectrum of the acceleration signal was 
averaged ten times to minimize noise effects. For the power spectrum, the analyzer can 
generate a second-order approximation for the pair of poles that correspond to the mode 
of vibration of the system. The poles are presented in the form 

s * D ±j f Hz. ( 41 ) 

From control theory, a pair of complex poles can be expressed as 

s » - 1 -C 2 Hz < 4 - 2) 


when 0 € ? £ 1 . 
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Comparing Equations (4.1) and (4.2), the natural frequency,/,, and damping ratio, f, 
can be solved from the analyzer output by 


/. « V/ 2+ J> : 


(4.3) 




(4.4) 


Notice that the analyzer output value / is the damped natural frequency. This value is 
actually of more use in this form when the modified command shaping technique is 
implemented in the controller. 

To obtain the desired frequency data, the workspace of RALF was divided up into 
ten degree joint increments and the power spectrum of the acceleration was taken at each 
location. Using Equations (4.3) and (4.4), the natural frequency and damping ratio of 
the first mode of vibration were calculated from the power spectrum output of the 
analyzer. 

4.3 Modal Analysis Results 

The fundamental damped natural frequencies calculated from the analyzer as a 
function of RALF’s joint coordinates are presented in Table 4.1. The fundamental 
natural frequencies of RALF were also calculated and are given in Table 4.2. Finally, 
the damping ratios were calculated from the analyzer output and are in Table 4.3. 
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Table 4.1 Fundamental Damped Natural Frequencies vs. Joint Coordinates 


fd (Hz) 

•2 (deg) | 



60 

70 

80 

90 

100 

107 


40 

5.17 

5.2 

5.17 

5.19 

5.46 

5.15 


50 

4.92 

4.99 

4.98 

4.87 

5.37 

5.24 


60 

4.75 

4.77 

4.84 

4.86 

4.99 

5.27 

0. 

(deg) 

70 

4.27 

4.26 

4.32 

4.28 

4.33 

wm 


3.79 

4.12 

4.11 

4.39 

4.37 

4.27 


90 

4.23 

4.29 

wm 

4.45 

4.51 

4.44 1 


100 

3.98 

4.07 

4.14 

4.25 

4.28 

4.33 


109 

wm 

3.8 

3.92 

3.97 

4.11 

4.07 | 


Table 4.2 Fundamental Natural Frequencies vs. Joint Coordinates 


|/.(Hz) 

(deg) ! 



60 

70 

80 

90 

100 

w 1 


40 

5.19 

5.21 

5.18 

5.21 

5.47 

u 

5.18 


50 

4.94 

5.01 

5.01 

4.9 

5.37 

5.24 


60 

4.76 

4.78 

4.85 

4.88 

5.01 

5.34 

0» 

(deg) 

70 

4.29 

4.29 

4.34 

4.30 

4.35 

4.42 

80 

3.8 

4.13 

4.12 

wm 

4.38 

4.28 


90 

4.25 

4.3 

4.41 

4.47 

4.52 

4.45 


100 

mm 

4.09 

4.16 

4.27 

UIH 

4.35 
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3.73 

3.83 

3.95 

3.99 

4.13 

4.12 


50 











































































































Table 4.3 Damping Ratios vs. Joint Coordinates 


I r 

d 2 (deg) 




60 

70 

80 

90 

100 

107 


40 

0.083 

0.071 

0.061 

0.085 

0.049 

0.107 


50 

0.084 

0.086 

0.116 

0.116 

0.026 

0.027 


60 

0.076 

0.076 

0.072 

0.079 

0.084 

0.166 

(deg) 

70 

0.105 

0.110 

0.093 

0.102 

0.091 

0.099 

80 

0.085 

0.083 

0.082 

0.076 

0.071 

0.065 


90 

0.093 

0.084 

0.078 

0.083 

0.062 

0.063 


100 

0.108 

0.098 

0.094 

0.086 

0.087 

0.092 


109* 

0.121 
P * 1 * 

0.129 

0.113 

0.109 

0.102 

0.149 


4.4 Curve Fitting of Experimental Data 

The damped natural frequency and damping ratio data is not useful for control 
purposes in tabular form. A look-up table to find needed values in a control algorithm 
can demand lots of precious computation time which slows down the sampling rate of the 
control system. Therefore, the data was curve fitted to shorten the computation time and 
permit a faster sampling rate. Appendix L discusses the least-squares curve fit of the 
experimental data. 
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The polynomial regressions that were derived for the damped natural frequency and 
damping ratio are 


//Ma> 55 - 


O.42120J + 1.8086 0^ - 2.1441 0 2 - 0.22290? 
+ 1.45400? - 3.63920J + 7.7830 


(4.5) 


C(0,,0 2 ) * - 0-42690? + 1.89180? -2.71590 2 - 0.05330? 
+ 0.34100? - 0.6148 0 t + 1.6762 . 


(4.6) 


To determine how well the equations approximated the experimental data, 
Equations (4.5) and (4.6) were evaluated at all of the joint positions. Each calculated 
value of the damped natural frequency, f c , was then divided by the experimental value 
of the damped natural frequency, f d , at the corresponding joint angles to create a 
normalized frequency, fjf d . The normalized frequencies as a function of joint 
coordinates is shown in Table 4.4. Using this normalized frequency, the correct number 
of impulses can be chosen based on the robustness criteria in Chapter IE. 

However, the damping ratio normalization generates the most error and is actually 
the determining factor for the number of impulses. Table 4.5 displays the normalized 
damping ratios, £/f, versus joint coordinates . This normalized data is very random in 
appearance which is a direct result of the way the analyzer linearly approximates the 
damping ratio. Using the largest normalized ratio of 4.5, Figure 3.13 suggests that the 
modified command shaping method should output four impulses to reduce vibration. 
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Table 4.4 Normalized Frequencies vs. Joint Coordinates 


fJU 

d 2 (deg) 


01 

(deg) 


60 

70 _ 

80 

90 

100 

107 

40 

0.992 

0.998 

1.02 

1.03 

0.990 

1.06 

50 

0.980 

0.980 : 

0.994 

1.03 

0.949 

0.978 

60 

0.961 

0.969 

0.970 

0.982 

0.970 

0.924 

70 

1.02 

1.04 

1.04 

1.07 

1.07 

1.06 

80 

1.10 

1.03 

1.05 

1.00 

1.02 

1.05 

90 

0.956 

0.956 

0.949 

0.956 

0.958 

0.978 

100 

0.988 

0.980 

0.981 

0.974 

0.982 

0.978 

109 

1.04 

1.03 

1.01 

1.02 

1.00 

1.02 


Table 4.5 Normalized Damping Ratios vs. Joint Coordinates 



0: (deg) S 



60 

70 

80 

90 

100 

107 


40 

1.63 

1.73 

2.12 

1.68 

0.30 

1.32 


50 

1.24 

1.06 

0.846 

0.959 

4.49 

4.13 


60 

1.12 

0.948 

1.10 

1.164 

1.16 

0.547 

0, 

fAatr\ 

70 

0.734 

0.582 

0.764 

0.825 

1.00 

0.831 

(deg; 

80 

0.917 

0.783 

0.870 

1.12 

1.28 

1.29 


90 

0.921 

0.865 

1.01 

1.11 

1.58 

1.45 


100 

0.913 

0.875 

0.980 

1.23 

1.29 

1.13 


109 

0.940 

0.781 

0.952 

1.11 

. 1.24 

0.800 


53 






























































































Now that the number of output impulses has been determined, the modified 
command shaping technique can be implemented to reduce vibration. The modified 
shaping method takes each input sample and replaces it with four output impulses that 
do not excite the first natural frequency of RALF. The required values of damped 
natural frequency and damping ratio needed to calculate the impulse amplitudes and their 
starting times are found from Equations (4.5) and (4.6). 
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CHAPTER V 


CONTROL IMPLEMENTATION AND RESULTS 


5.1 Control Structure 

The original input shaping method devised by Singer was entirely feedforward in 
design [38]. The desired (x,y) coordinate positions were transformed into desired joint 
angles using a Jacobian before applying the input shaping algorithm. This feedforward 
scheme provides little robustness to noise disturbances or to model uncertainty. To 
overcome these problems, Hillsley and Yurkovich [18] applied a composite control 
strategy which utilizes input shaping with a feedback scheme. 

The control design for RALF uses an adaptive, proportional plus derivative (P.D.) 
feedback strategy derived by Yuan [49] with the addition of modified command shaping. 
Yuan’s control algorithm was chosen because it compensates for unmodeled modes and 
nonlinearities of the system. The modified command shaping technique is designed to 
eliminate the first natural frequency of RALF discussed in Chapter IV. The block 
diagram of the comprehensive control system is shown if Figure 5.1. Since the input to 
the control system is desired angles, it is easily implemented with the master arm which 
allows the user to specify a desired end-position while it records desired joint angles. 
This control implementation also allows other pre-computed trajectories to be input into 
the system for means of comparison. 
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Figure 5.1 Block Diagram of Control System 


5.2 Input Trajectories 

To demonstrate the reduction in tip vibration, the modified command shaping 
technique is compared to the original P.D. controller for different input trajectories. For 
these comparisons, the trajectories are pre-computed to ensure an equivalent basis for 
comparing vibration reduction. The principal trajectory is a three-foot diameter circle 
located above the first axis with- a completion time of nine seconds. Figure 5.2 shows 
the location of the circle relative to the workspace of RALF. This location was chosen 
because there is a large enough variation in system natural frequency to change the 
discrete-time period, deln. A sinusoidal perturbation signal with variable frequency is 
then added to the radius parameter of the circle to induce vibration into the system. 

The vibration of the tip is recorded using two different methods. The first 
method is very similar to the way the frequency data was obtained in Chapter IV. The 
robot is commanded to follow the desired circle trajectory eight times which eliminates 
noise effects and allows reliable averaging of the data. The analyzer records the 
transverse acceleration response at the tip and computes the power spectrum of the data. 
The magnitude of the frequency response is the root-mean-square of the acceleration and 
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Figure 5.2 Principal Circle Trajectory 


is displayed in units of decibels (dB). Recall the definition of decibels is 

number of decibels - 10 log 10 A (5.1) 

where A is the amplitude in question. The frequency response for each control strategy 
is then compared to determine which controller has the smaller amplitude of vibration. 

The second method is a visual approach to examine the tip vibration of the robot. 
A light-emitting diode (LED) is attached to the end-point of RALF and then one 
revolution of the desired circle trajectory is commanded. By leaving the aperture on a 
35 mm camera open, the actual end-point path is recorded. Of course, special measures 
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are made to ensure that the LED is the only light source imprinted on the negative. At 
the end of the trajectory, the flash is triggered to reveal the robot configuration. 

5.2.1 Trajectory One - Circle With No Perturbation 

The first trajectory is the principal circle with no perturbation added. Figure 5.3 
shows a comparison of the frequency response between the P.D. control and Singer’s 
original input shaping substituted for the modified command shaping. The magnitude of 
the frequency response for the input shaping method is greater than that for the P.D. 
routine. 



Figure 5.3 Frequency Response of RALF 
P.D. vs. Input Shaping 
No Perturbation 


58 



The increase in magnitude of the frequency response for the input shaping method 
at the 4.8 Hz frequency value denotes a vibration of the system at the first natural 
frequency. This vibration is provoked by the inability of the input shaping to handle 
varying system parameters. The sharp peak in the frequency response at the 10 Hz 
frequency value is the second natural frequency. However, the modified command 
shaping technique is only designed to reduce the vibration at the first natural frequency. 

The induced vibration caused by the input shaping method can be seen visually 
in Figure 5.6. The path followed by the end-point using the ordinary P.D. controller is 
shown in Figure 5.5. The actual path followed by the end-point of RALF is not a 
precise circle because the P.D. controller does not drive the steady state joint error to 
zero. 

By implementing the modified command shaping technique, the induced vibration 
is eliminated. This result is evident in the frequency response of Figure 5.4 and the 
picture given in Figure 5.7. From Figure 5.4, the magnitude is reduced by 20 dB at the 
system natural frequency of 4.8 Hz. This results in a vibration that is only 1 % of the 
amplitude of the original P.D. vibration for this specific frequency. Comparing Figure 
5.5 and Figure 5.7, the modified command shaping technique generates almost the same 
tip motion as the P.D. controller even though the 4.8 Hz frequency component is 
missing. 
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Figure 5.4 Frequency Response of RALF 

P.D. vs. Modified Command Shaping 
No Perturbation 



Figure 5.5 Picture of Tip Motion Using P.D. Control 
No Perturbation 
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Figure 5.6 Picture of Tip Motion Using Input Shaping 
No Perturbation 



Figure 5.7 Picture of Tip Motion Using Modified 
Command Shaping - No Perturbation 
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5.2.2 Trajectory Two - Circle With 1 Hz Perturbation 

This trajectory, shown in Figure 5.8, is the principal circle with a 1 Hz sine wave 
with an amplitude of 1.5” riding on the radial component of the circle. This trajectory 


0— Trajectory 
** **** ** Ht 



Figure 5.8 Circle Trajectory with 1 Hz 
Perturbation 


should contain nine "bumps* around the circle since the period is nine seconds. The 
frequency response of the modified command shaping versus the P.D. control is shown 
in Figure 5.9. Since the command shaping technique was not designed to eliminate 1 Hz 
vibration, the two control schemes show comparable results for this frequency range. 
However, the modified command shaping reduced the magnitude of vibration by 18 dB 
at the system natural frequency value of 4.8 Hz. This results in a vibration that is 1.6% 
of the amplitude of the original P.D. vibration at this particular frequency. 
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Figure 5.9 Frequency Response of RALF 

P.D. vs. Modified Command Shaping 
1 Hz Perturbation 


5.2.3 Trajectory Three - Circle With 4.8 Hz Perturbation 

The purpose of this trajectory, given in Figure 5.10, is to excite the first natural 
frequency of the system. Figure 5.11 displays the frequency response comparison for 
this trajectory. The difference in magnitude is 32 dB at the system natural frequency of 
4.8 Hz which corresponds to 0.06% of the original P.D. vibration amplitude for this 
particular frequency. The visual effects are even more impressive. Figure 5.12 shows 
the tip motion for the P.D. control effort while Figure 5.13 displays the tip motion for 
the modified command shaping. At least visually, the modified command shaping 
method appears to completely eliminate the vibration. 
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Figure 5.10 Circle Trajectory with 4.8 Hz 
Perturbation 



Figure 5.11 Frequency Response of RALF 


P.D. vs. Modified Command Shaping 
4.8 Hz Perturbation 
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Figure 5.12 Picture of Tip Motion Using P.D. Control 
4.8 Hz Perturbation 



Figure 5.13 Picture of Tip Motion Using Modified 
Command Shaping - 4.8 Hz Perturbation 
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5.2.4 Trajectory Four - Circle with 10 Hz Perturbation 

The final circle trajectory, displayed in Figure 5.14, demonstrates the ability of 
the modified method to reduce system frequency vibration for input signals that contain 
higher frequency components. Figure 5.15 shows a reduction in magnitude at the 


Desired Tropctory 
Mew i w fne<«*y i ojohz 



Figure 5.14 Circle Trajectory with 10 Hz 
Perturbation 


4.8 Hz frequency location of 9 dB. This results in a vibration that is 12.6% of the 
amplitude of the original P.D. vibration for this particular frequency. 


66 




Figure 5.15 Frequency Response of RALF 

P.D. vs. Modified Command Shaping 
10 Hz Perturbation 


5.2.5 Trajectory Five - Pseudo-Step 

A step input is often given as a test input because it theoretically contains all the 
frequency components. To simulate a step, a cycloidal motion in joint space was created 
with a duration time of 0.2 seconds [30]. Using this pseudo-step input, a ten degree step 
was simultaneously given to each joint and a time record was taken of the transverse 
acceleration response at the tip. Figure 5.16 shows the desired tip motion with the 
starting position at (-4.7,142.9) and Figure 5.17 displays the transverse acceleration 
response. Clearly, the amplitude of the time response for the pseudo-step input is 
reduced using the modified method. However, notice the delay in the modified command 
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Figure 5.16 Pseudo-Step Trajectory 



Figure 5.17 Pseudo-Step Response of RALF 

P.D. vs. Modified Command Shaping 
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shaping method’s pseudo-step response. This delay is present because the pseudo-step 
trajectory holds the robot in the starting joint configuration for ten clock cycles. Since 
the sampling rate is 50 Hz, this time delay is 0.2 seconds. Although the robot is 
supposed to remain stationary during the hold time, a slight error signal is developed in 
the control algorithm which causes the robot to move. This tip motion is large enough 
for the amplitude of the transverse acceleration to trigger the analyzer. This slight tip 
motion is verified in the Figure 5.18. When the modified command shaping method is 
implemented, the amplitude of the small error signal is reduced so that the transverse 
acceleration is not large enough to trigger the analyzer. Therefore, the delay experienced 
by the modified command shaping method is actually the hold time of the pseudo-step 
trajectory. 

Figures 5.18 and 5.19 show the pseudo-step responses for each of the control 
schemes. The slight tip motion during the hold portion of the trajectory is very evident 
and is even more pronounced in the modified command shaping technique case. The 
reduction in vibration of the modified command shaping method is apparent from the 
smoothness of the trajectory. However, the tip position overshoots the desired final 
position. The feedback nature of the control structure together with an underdamped 
system is believed to generate this undesirable overshoot. Had the method been 
implemented in a feedforward manner as Singer did, this overshoot might be eliminated. 
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Figure 5.18 Picture of Tip Motion Using P.D. Control 
Pseudo-Step Input 



Figure 5.19 Picture of Tip Motion Using Modified 
Command Shaping - Pseudo-Step Input 
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CHAPTER VI 


CONCLUSIONS 


6.1 Summary 

Singer’s original input shaping idea was introduced and explained using the 
equation of motion for a simple second-order system. The constraint equations needed 
to yield a vibrationless system were discussed and then multiple impulse amplitudes were 
derived for several constraint conditions. A vibration error expression was derived based 
on the ratio of the multiple impulse response to the impulse response of a second-order 
system. The residual vibration for each case was graphed to determine the effects of 
deviation in system parameters. The input shaping method was shown to actually 
increase the residual vibration in certain situations when the system parameters are not 
accurately specified. 

The main intention of this research was to reduce end-point vibration in a 
teleoperated system while preserving the usefulness of the system motion. Results 
verified that the input shaping method can actually induce vibration in systems that have 
varying parameters. Therefore, a modified command shaping technique was developed 
to alleviate this problem. By varying input perturbations to trajectories in a range of 
possible user inputs, the modified command shaping technique proved to reduce vibration 
at the system’s first natural frequency. The reduction in amplitude varied from 87.4% 
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for the high frequency perturbation to 99.94% for a perturbation frequency near the 
system resonance. The acceleration response of the command shaping to a pseudo-step 
input was smaller in amplitude than the regular P.D. control. However, the pseudo-step 
response of command shaping displayed overshoot which may be undesirable in many 
end-point positioning tasks. 

6.2 Contributions 

The major contributions of this research are the implementation of teleoperation 
on RALF, the new perspective of residual vibration, the visualization of tip motion and 
the derivation of a modified command shaping technique. 

The implementation of teleoperation on RALF allows users to specify any end- 
point trajectory within the workspace. The ability to perform teleoperated experiments 
such as cutting and end-point bracing are now possible. 

The new definition for residual vibration gives a clearer representation of the 
input shaping method. It demonstrates that the input shaping method may not be 
effective in reducing residual vibration in all cases. 

The visualization method for viewing tip motion using a LED permits easy 
evaluation of tip vibration. The method can also estimate the ability of a system to 
follow desired tip trajectories. 

Finally, the modified command shaping technique expands the use of the original 
input shaping idea to systems with varying parameters. The induced vibration caused 
when the input shaping method is applied to a variable frequency system is eliminated 
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by the command shaping method. By realizing a change in the discrete-time period, the 
impulse input is modified to eliminate the possibility of producing the vibration. 

6.3 Future Work 

The main concern to be addressed in future research is the elimination of 
overshoot in the step response of RALF. To have an effective teleoperated system, the 
end-point must be positioned without overshoot. This could be accomplished using the 
modified command shaping technique in a feedforward arrangement that is combined 
with a feedback control scheme. Another possibility is the implementation of inverse 
dynamics with the modified method to accurately position the end-point of RALF. 

The second goal is to implement the modified method to eliminate the second 
natural frequency of RALF. This would require a fester control system than the 
MicroVAX can provide. By using a PC, faster computation and A/D rates are possible 
and the elimination of the first two natural frequencies could be achieved. 

Finally, the effects of the modified command shaping technique on the stability 
of the overall control structure will be investigated. The time delays produced by this 
shaping method must be better understood to develop a control scheme that has accurate 
end-point tracking capabilities without overshoot. 
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APPENDIX A 


RIGID BODY TRANSFORMATIONS FOR RALF 

Forward kinematics of rigid bodies is used to determine the position and 
orientation of the end-effector on a manipulator. Inverse kinematics, on the other hand, 
with the problem of rinding all the possible joint configurations given the end- 
effector position and orientation. Questions often arise with the uniqueness of a given 
solution, or even if one exists, with the inverse case. 

By using matrix transformations, kinematic equations relating the end-point of a 
manipulator to its base can be derived. The process involves representing each degree 
of freedom of a manipulator with an individual coordinate system. Each matrix 
transformation relates one coordinate system to the previous one. The overall matrix 
transformation from the end-point to the base of a robot is obtained by multiplying all 
of the matrix transformations in the proper order. Most robotic texts give a thorough 
discussion of matrix transformations [11,44]. 

Although RALF is not a rigid body, the forward and inverse kinematics can be 
derived to give an approximate solution for the end-point location and orientation or the 
joint configurations. Consider Figure A.l which shows the rigid body coordinate 
systems used for RALF. 
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1 t 
2 1 


(A.2) 


cos(0 2 ) -sintSj) 0 L x 
sin(0 2 ) cos(0 2 ) 0 Ly 

0 0 10 
0 0 0 1 



0 0 1 , 

1 0 0 

0 1 0 

0 0 1 


(A.3) 


The overall matrix transformation from the end-point to the base of RALF is found by 
multiplying the individual transformations as follows 

Jr* Jr Jr Jr. < A - 4 > 


The resulting overall matrix transformation is 



cos(0,+0 2 ) 

sin(0,+0 2 ) 

0 

0 


-sm(0j-*-0 2 ) 

cos(0,+0 2 ) 

0 

0 


0 I^co^Oj+OjJ+LjCosCO^-IjSinCOj) 

0 IjSinCOj+Oj) ♦IjanCOj) +L 3 cos(B l ) 

1 0 

0 1 


(A.5) 
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Recall that the fourth column of the overall matrix transformation is the position vector. 
Therefore, the (x,y) coordinates of the end-point are given in base coordinates by 

x * LjCosfSj +0 2 ) + Lj cos(0j) - I 3 sm(0,) (A. 6 ) 

y - LjSmCQj +0 2 ) + I, sinCOj + 1,008(0,) . (A.7) 

Given the two joint angles, the (x,y) coordinates of the aid-point can be calculated using 
Equations (A. 6 ) and (A.7). 

Now consider the inverse problem. Taking the sum of the squares of Equations 
(A. 6 ) and (A.7) yields 

x 2 * y 2 - i} + 1% * 1% * 2L,I 2 cos(0 2 ) + ILzLsSmiBJ . ( A -8) 

Define K by 

< A - 9 > 

so that Equation (A. 8 ) becomes 

K « 21 , 12008 ( 62 ) + 2L 2 I 3 sin(6 2 ) . (A. 10) 


77 



Solving for cos(0j) from Equation (A. 10) yields 


cos(0 a ) 


!T-2I^LjSin(0 3 ) 

2l^ 


(A. 11) 


Notice that Equation (A. 11) can not be explicitly solved for 0 2 . But recall the 
trigonometric identity 


sin 2 (0 2 ) + cos 2 (0 J ) - 1 . 


(A. 12) 


Squaring Equation (A. 11) 


co^COj) = 


¥} -4^Ljl3sin(0 2 ) +4I?l£sm 2 (0 2 ) 


(A. 13) 


and substituting into Equation (A. 12) yields 


sin 2 (0 2 ) = 


4I?lj & +AKL 2 I,sin(0 2 ) -4ljl?sm 2 (0 2 ) 

Al^li 


(A. 14) 


Equation (A. 14) can be simplified to 


sin(0 2 ) 


ITI^±Z^4lj(i;+lj)-Jg 2 

21,(1? *li) 


(A. 15) 
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Since d 2 can not be less than zero, Equation (A. 15) becomes 


21,(1? +4) 


Substituting Equation (A. 16) into Equation (A. 11) produces 


cos(0 2 ) = 1 — 

22 , 0 ? * 1 ?) 


Using Equations (A. 16) and (A. 17), the solution for d 2 becomes 


0 2 * atan2 


( KLj * L t \j + 1?) -2T 2 ] 




To solve for Equations (A.6) and (A.7) are expanded to produce 
x * (L 2 cos(0 2 )+L 1 )cos(0 1 )-(I 2 sin(0 2 ) +L 3 )sm(Q 1 ) 


y - (^Sm^+^COS^+^COS^+I^Sintf,). 


(A. 16) 


(A. 17) 


(A. 18) 


(A. 19) 


(A. 20) 
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Equations (A. 19) and (A.20) are solved for sin(0,) and cos(0 t ), which are 


-xfl^sm^) •*•!,) +y(I,cos(8 2 ) +1,) 
lZ+l}+lZ +2L 2 (L l cos(d 2 ) +LjSin(0j)) 


008(0!) - 


xiLjCCsQi) +!,) *y^2 ,h ^a) ihl 
lZ+li+lZ+ 21,(1, cos( 0 2 )+l 38 m( 0 2 )) 


Using Equations (A.21) and (A.22), the solution for 0, becomes 


6i 


atan2 


f -jc(I,sin(0 2 ) +LJ -«-y(I a cos(0 2 ) •*•!,) \ 

^ x(I 2 cos(0 2 )+L,)+y(I 2 sm(0 2 )+I 3 ) J’ 


(A.21) 


(A.22) 


(A.23) 
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So, knowing the (x,y) position of the end-point of RALF, the joint coordinates can be 
computed using the following three equations: 

K-* -il-li-lZ (A - 24) 


e, 


■ atari! 


< KL r L 3 {7iiU^y^ j 


(A.25) 




atan2 


' -x(LjSm(0 2 ) +L ,) +y(l5cos(0 2 ) +L X ) \ 
k x(L 2 cos(0 2 )+L 1 )+y(L 2 sm(Q 1 )+L 3 ) j 


(A.26) 
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APPENDIX B 


DIMENSIONS OF MASTER ARM 


Figure B. 1 shows the dimensions of the master arm used by the operator to input 
desired joint commands. Potentiometers, listed in Appendix C, are placed at each joint 
of the master arm to translate the joint commands input by the user into actual joint 
motions of RALF. 
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Figure B. 1 Dimensions of Master Arm 
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APPENDIX C 


EQUIPMENT LIST 


ELECTRONIC COMPONENTS 


Device: MicroVax II 
Model No.: VS21W-A2 
Serial No.: WF61305805 
Company: Digital Equipment Co. 

Device: Real-Time Clock 
Model No.: DT2769 
Serial No.: 187824-C453 
Company: Data Translation 

Device: Analog I/O System 
Model No.: DT2785 
Serial No.: 188872-C496 
Company: Data Translation 

Device: Dynamic Signal Analyzer 
Model No.: HP3562A 
Serial No.: 2502A00718 
Company: Hewlett-Packard Co. 

Device: Disk Storage Unit 
Model No.: HP9122D 
Serial No.: 2518A44227 
Company: Hewlett-Packard Co. 

Device: Supply/ Amplifier - Joint one of RALF 
Model No.: BOP36-1.5M 
Serial No.: F79808 
Company: KEPCO 
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Device: Supply/Amplifier - Joint two of RALF 
Model No.: BOP36-5M 
Serial No.: F105062 
Company: KEPCO 

Device: Potentiometer - Joint one of Master 
Model No.: 6637 
Company: Bourns 

Device: Potentiometer - Joint two of Master 
Model No.: 6637 
Company: Bourns 


SENSORY COMPONENTS 


Device: LDT Position Sensing System - Joint one of RALF 
Model No.: 011020050100 
Serial No.: 20658-02-002P 
Company: MTS Systems Corp. 

Device: LDT Position Sensing System - Joint two of RALF 
Model No.: 01 1020050100 
Serial No.: 20658-02-001P 
Company: MTS Systems Corp. 

Device: Analog Output Module - Joint one of RALF 
Model No.: 0110200503105001 
Serial No.: 20658-05-001P 
Company: MTS Systems Corp. 

Device: Analog Output Module - Joint two of RALF 
Model No.: 0110200503105001 
Serial No.: 20658-04-001P 
Company: MTS Systems Corp. 

Device: ICP Accelerometer 
Model No.: 308B 
Serial No.: 10430 
Company: PCB Piezotronics 
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Device: Power Unit 
Model No.: 480D06 
Serial No.: 5163 
Company: PCB Piezotronics 


HYDRAULIC COMPONENTS 


Device: A.C. Motor 
Model No.: 2U2100 
Serial No.: D-75 
Company: Delco 

Device: Vickers Variable Volume Piston Pump 
Model No.: F3-PVB20-FRS-20-C-11 
Company: Sperry Rand Corp. 

Device: Hydraulic Valve - Joint one of RALF 
Model No.: 73-102A 
Serial No.: 144 
Company: Moog, Inc. 

Device: Hydraulic Valve - Joint two of RALF 
Model No.: 73-102A 
Serial No.: 147 
Company: Moog, Inc. 

Device: Hydraulic Cylinder - Joint one of RALF 
Model No.: H-PB-2 
Serial No.: 3778 1-J 
Company: Atlas Cylinder Corp. 

Device: Hydraulic Cylinder - Joint two of RALF 
Model No.: N2C-3.25x40 
Serial No.: 5C8205-065-1B 
Company: Hydro-Line Mfg. Co. 
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APPENDIX D 


IMPULSE RESPONSE OF A SECOND-ORDER SYSTEM 



F(t) 


Figure D.l Spring-Mass-Damper System 


Recall that Newton's second law is 

Y F=ma . 


/ 


(D.l) 
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Newton’s second law applied to the simple system shown in Figure D.l yields 


which can be rearranged to get 

d}x dx . 

m — — + c — + kx - F(t ) . 

dx 1 & 


Dividing through by the mass, the system equation becomes 

dh + S_dx + _ F(t) 

dt 1 m dt m m 


By definition, the natural frequency of a second-order system is 



and the damping ratio is 



(D.2) 


(D.3) 


(D-4) 


(D.5) 


(D.6) 
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Substituting into the system equation, 


— + 2C“.^ * ". J * * —TO. 

dt c 


(D.7) 


By defining an impulse force as 


fit) = 


F(t) 


(D.8) 


the system equation becomes 


^ ♦ «. J * - ®„ 7 <») 


<fr 




(D.9) 


Let the impulse force have the form 

f(t)=A6<t-t 0 ) (DIO) 

where A is an amplitude and 5(f) is the Dirac delta function. A Dirac delta function is 
defined mathematically in Equation (D.ll) and is displayed in Figure D.2 on the 
following page. 

J^6(f-f 0 )«l (D- 11 ) 
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Figure D.2 Dirac Delta Function 


From Figure D.2, it is apparent that the Dirac delta function has an amplitude of 
l/(time). The impulse function f(t), given in Equation (D.10), must have units of 
length. This results in the amplitude A having units of (length *time). It is sometimes 
very helpful to conduct unit analysis of this type to better understand the physics of the 
problem. 

The new system equation becomes 

^ - 2(u — ♦ o’* - a.’A »(r (D.12) 

dp dt 

Recall the definition of the Laplace Transform 

F(s) - J“/(r)e'"dr. (D* 13 ) 
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Taking the Laplace Transform of the system equation assuming zero initial conditions 
yields 


s 2 X(s) * 2 Cu m sX(s) «■ u 2 X(s) - Au.V " 0 


(D.14) 


which can be rearranged to get 


X(s) = 


s 2 + + «* 2 


(D.15) 


Partial fraction expansion of the equation yields 


X(s) = 


JAu e’* ,#/ 




s + Cw, + ju H fi ~ C 2 s + («„ - C 2 J 


(D.16) 


The inverse Laplace Transform of this expression is 

*( f ) = (o >17 ) 

2/W 2 


91 



which is just 


yw 5 ^ v 


Recalling the Euler representation for sine is 


sin(o) - 




the second-order system response to an impulse force becomes 


x(r) « — " sanlu^Jl - f 2 (* -f 0 )) 

/W 5 
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APPENDIX E 


SUMMATION OF TWO SINUSOIDS 


The following two sinusoidal responses, 

XjCO-^sinCaf+^i) (E*l) 

Xj (r) * £ 2 sin(a t + 4> 2 ) (E*2) 

can be linearly combined to form a total response, x(r), where 

x(r)=B 1 sin(af + <J> 1 )+£ 2 sin(af+4>2) • (^.3) 

Recall that 

sin(a+l>) !B sin(a)cos(b) +sin(b)cos(a) . 

So the total response becomes 

x(r) =B 1 [sin(at)cos(<J) 1 ) +sin(4> 1 )cos(af)] +5 2 [sin(ar)cos(4> 2 ) +sin(<t> 2 )cos(af)] (E.5) 
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which is just 


X (f) - [Bjcos^) + R 2 cos(<t> 2 )] sin(ar) + [B 1 sin(<t> 1 ) +B 2 sin(<t> 2 )]cos(af) . 


Recall that the Euler representations for sine and cosine are 

Ju _ M -Ju . - -ju 

sin(w) (E.7) cos(o) * . 

2 / 2 


Substituting, 


x(t) = [J5 1 cos(4> 1 ) +fl 2 cos(4> 2 )] 




V 


+ [BjSinWj) +B 2 sin(<|> 2 )] 




Rearranging the coefficients of the exponential toms. 


x(f) * [(fi,cos(4>,) +B 2 cos(4» 2 )) +y(fl,sin(4> 1 ) +R 2 sin(<|> 2 ))]^ r 


- [(Bjcostoj) +B 2 cos( 4 > 2 )) -/(RjSin^j) +fl 2 sin(<|> 2 ))]4^- • 

A J 


Recall that any complex number can be written in the form 


z=re J * 


where r is the magnitude and 6 is the phase. 
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(E.6) 


(E.8) 


(E.9) 


(E.10) 


(E.11) 


■ 


Using this idea, x(i) can be expanded as follows 

’ Jtm-j *i**+i>**i aiDi *> \] . f 

x(t) = ^/[B 1 cos(4> 1 ) +B 2 cos(4> 2 )] 2 + [BjSinC^j) +B 2 sinOfr 2 )] 5 e '* l "^ 4|> — 

\a,oc*(* l )*a^(^/ e 7 
J 2 j 

(E* 12) 


and simplifying, 

x(f) - [ ^[BjCos(<|>,) +B 2 cos(4> 2 )] 2 +[B 1 sin(4> l ) +£ 2 sin($ 2 )] 5 ] 

(E.13) 
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This result can be rewritten as 


x(/) ■ g^sintef + f ) 


(£•15) 


where 


g^- ^[gjCosC^) +g 2 cos(<|> 2 )] 2 + [gjSm(4>,) +g2sm(4> 2 )] 2 


(E.16) 


g x cos(4»,) +g 2 cos(<t> 2 ) ; 


(£•17) 
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This result can now be generalized. The total response, x(f), which is the linear 
combination of N sinusoidal responses can be written as 

x(jt) • B^smiat + +) (E- 1 8) 


where 


B 


amp 


N 


E B k cos(<fr t ) 




(E.19) 


* 


/ N \ 

Y B k sm(4> t ) 


tan 


-i 


*■1 


[E^cos^J 


(E.20) 
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APPENDIX F 


CALCULATIONS FOR TWO IMPULSES 


The constraint equations that must be satisfied are 

i^cosMj) +£ 2 cos(4> 2 ) » 0 


^sinWi) + * a sia(<M ■ 0 


where 


B k 


g -C «.(»-*«) 

fTT 1 


and 


4> t » - * 


(F.l) 


(F.2) 


(F.3) 


(F.4) 


Notice that there are four unknowns ( A„ A 2 , t 01 and and only two constraint 
equations. The other two constraints are the starting time of the first impulse and the 
normalization of the amplitudes. 


Pick t 01 = 0 so that Equations (F. 1) and (F.2) become 


Bj + B 2 cos(<|> 2 ) - 0 


(F.5) 


B 2 sin($ 2 ) » 0. 


(F.6) 


The two solutions that exist for Equation (F.6) are either B 2 - 0 or 

<|> 2 » ntc 


(F.7) 


where n = ± 0,1, 2, 3,... . Avoiding the trivial solution, pick n = - 1 so that 


♦2 


- it . 


(F.8) 


From Equation (F.4), the second impulse occurs in time at 


t 


02 



(F.9) 


Substituting Equation (F.8) into Equation (F.5) yields 




(F.10) 
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At t - t m , Equation (F. 10) produces 


i ^l fa) ii f C ~c«». (<<!-<«> 

JTT 1 v'w 5 


or 


Aj « A x t 


c« 

/T 7 ? _ 


If the amplitudes are normalized so that 


2 


5>, 


i 


and let 


M * e 


c* 


Vl-C* 

» 


(F.U) 


CF.12) 


(F.13) 


(F.14) 
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then the amplitudes of the two impulses and the times at which they occur are 




1 

1 + Af 


(F.15) 


f 01 - 0 (F.16) 


A , - (F IT) 

2 1 + Af 


*02 



(F.18) 


Now that the amplitudes and times of the impulses are known, the theoretical impulse 
response of a second-order system can be found using the results from Appendix D. 
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APPENDIX G 


CALCULATIONS FOR THREE IMPULSES 


The constraint equations that must be satisfied are 

+£ 2 cos(4> 2 ) +fl 3 cos(4> 3 ) ■ 0 (G.l) 

B,sm(4»,) + £ 2 sm(4> 2 ) +B 3 sm($ 3 ) ■ 0 (G.2) 

cos (4^) * fl 2 f 02 cos(<|> 2 ) +fi 3 f 03 cos(<)> 3 ) - 0 (G.3) 

+ 5 2 r 02 sin(4> 2 ) + £ 3 f 03 sm(4> 3 ) = 0 (G.4) 


where 

B - A * n * c -C«h«-W 

* /w* 
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and 


♦» • -<Vo*Vi -c J - 


(G.6) 


Notice that there are six unknowns (A lt A 2 , A Jt t 01 , t m and t &) and only four constraint 
equations. The other two constraints are the starting time of the first impulse and the 
normalization of the amplitudes. 

Pick t 01 = 0 so that Equations (G.l) through (G.4) become 

Bj + B 2 cos(4> 2 ) + B 3 cos(4> 3 ) - 0 (G- 7 ) 

B 2 sm(4> 2 ) +B 3 sin(<t> 3 ) - 0 (G.8) 

B 2 f 02 cos(4> 2 ) +B 3 r 03 cos(<|> 3 ) * 0 (G.9) 


B 2 r 02 sin(4> 2 ) +B 3 r 03 siu(<|> 3 ) - 0. 


(G.10) 


From Equation (G.8), 


B 2 sin(<J> 2 ) - -B 3 sin(4> 3 ) 


(G.ll) 
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which is substituted into Equation (G. 10) to produce 




The three solutions that exist for Equation (G. 12) are 


*3 


- 0 


*03 ~ *02 


- 0 


4> 3 ■ nit 


where n = ± 0,1, 2, 3,... . Avoiding the trivial solutions, pick n = - 2 so that 

4*3 = -2 It. 


From Equation (G.6), the third impulse occurs in time at 


2 it 


*03 


W .v^l-C 2 


(G. 12) 


(G.13) 


(G.14) 


(G.15) 


(G.16) 


(G.17) 
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Substituting Equation (G.16) into Equation (G.8) yields 


B 2 sin(4> 2 ) - 0. 


(G.18) 


The two solutions that exist for Equation (G.18) are 



(G.19) 


<J> 2 - nit 


(G.20) 


where n - ± 0,1,2, 3,... . Avoiding the trivial solution, pick n = -1 so that 



(G.21) 


From Equation (G.6), the second impulse occurs in time at 


*02 = 


It 


1 ' C 2 


(G.22) 


Now that the times at which the impulses occur is known, the amplitudes of the three 
impulses can be found. Substituting Equations (G.16) and (G.21) into Equations (G.7) 
and (G.9) produces 

b,-b 2 + b 3 * o ( G - 23 ) 
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(G.24) 


-B 2 t 02 + B 3 t 03 -0 


From Equations (G.17) and (G.22), the relationship between the times is 


*03 * 2*02 * 


Substituting this time relationship into Equation (G.24) yields 


* 2-2 *1 


*3 -*!• 


At r = Equation (G.26) produces 


***** c -Cu.d w -i M ) x 




/I 7 ? 


or 


A 2 - 2 A x t 


(G.25) 


(G.26) 


(G.27) 


(G.28) 


(G.29) 
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At t = Equation (G.27) produces 


T 1 


or 


\ 3 


i4j« 


-== 
V'l-f* . 


If the amplitudes are normalized so that 

E4 - i 

i-1 


and let 


M * e 


c« 


Vi-c 2 

t 


(G.30) 


(G.31) 


(G.32) 


(G.33) 
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then the amplitudes of the three impulses and the times at which they occur are 


A. - (G.34) r « 0 (G.35) 

1 1+2M + M 2 




2Af 

1 + 2M + M 2 


(G.36) 


*02 



(G.37) 


A 3 


1 + 2M + M 2 


(G.38) 


*03 


2 ic 

1 ~ c 2 


(G.39) 


Now that the amplitudes and times of the impulses are known, the theoretical impulse 
response of a second-order system can be found using the results from Appendix D. 
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APPENDIX H 


CALCULATIONS FOR FOUR IMPULSES 


The constraint equations that must be satisfied are 

BjCosC^j) + B 2 cos(<J> 2 ) + B 3 cos(4> 3 ) + £ 4 cos(<|> 4 ) * 0 (H.l) 

5^(4),) + B 2 sin(4> 2 ) + B 3 sm(<|) 3 ) + B 4 sin(<|> 4 ) - 0 (H.2) 

B x t 0l cosC^) + B 2 f 02 cos(4» 2 ) + B 3 t 03 cos(4> 3 ) +B 4 t 0A cos(*J ■= 0 (H.3) 

B t f 0 jSinC^j) + B 2 f 02 sin(4> 2 ) + B 3 t oi sin(<|> 3 ) + B 4 f 04 sin(<t> 4 ) * 0 (H.4) 

Bj^cosC^,) ♦ B 2 r 0 2 2 cos(4> 2 ) ♦ B 3 fJ 3 cos(4> 3 ) «■ B 4 & cos(4> 4 ) - 0 (H- 5 ) 

B^sinC*,) ♦ B 2 f 0 2 2 sin(4) 2 ) ♦ B 3 r 2 3 sin(4> 3 ) ♦ B 4 t 2 04 sin(4> 4 ) - 0 («.6) 
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where 


B k 


/w 5 


(H.7) 




Notice that there are eight unknowns ( A l% A 2 , A Jf A 4 , t 0I , t m , and /<*) and only six 
constraint equations. The other two constraints are the starting time of the first impulse 


and the normalization of the amplitudes. 

Pick t 0l = 0 so that Equations (H. 1) through (H.6) become 

+ * 2 cos(4> 2 ) + £ 3 cos(4> 3 ) + B 4 co$(4> 4 ) - 0 (H.9) 

B 2 sin(4> 2 ) + B 3 sin(4> 3 ) + U 4 sin(<|> 4 ) « 0 (H.10) 

B 2 r 02 cos(4> 2 ) + B 3 f 03 cos(4> 3 ) + £ 4 r M cos(+ 4 ) » 0 (H.ll) 

£ 2 f 02 sin(<j> 2 ) + £ 3 f 03 sin(4> 3 ) + f 4 r M riit(4 4 ) * 0 (H. 12) 


110 


£ 2 fJ 2 cos(<|> 2 ) + £ 3 fJ 3 cos(4> 3 ) + S 4 fJ 4 cos(<j> 4 ) - 0 (H- 13 ) 

£ 2 *J 2 sin(4> 2 ) + £ 3 fJ 3 sm(4> 3 ) + B^sani * 4 ) * 0 • (H,14) 

From Equation (H.10), 

£ 2 sin(4> 2 ) - - [£ 3 sin(<j> 3 ) + £ 4 sin(<|> 4 )] (H* 15 ) 

which is substituted into Equations (H.12) and (H.14) to produce 

* 0 <H16) 

fi 3 (( 0 J J -4)sm(*,)4 B 4 «M-*2 2 )“>(+4) - 0- (H17) 

Solving for B, from Equation (H. 16) and substituting into Equation (H.17) yields 

B 4 (r 0 , - r„)sin(* 4 ) - 0. (H.18) 

The three solutions that exist for Equation (H.18) are 

B 4 - 0 (H.19) 
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(H.20) 


*04 ~ *03 " ® 


<(> 4 ■ nit 


where n * ± 0,1, 2,3,... . Avoiding the trivial solutions, pick n * - 3 so that 

<|>4 - - 3 ti . 


From Equation (H.8), the fourth impulse occurs in time at 


t 


04 


3 « 


Substituting Equation (H.22) into Equation (H.16) yields 

B 3 (t 03 -t 02 ) sin ( < J > 3) * °- 


The three solutions that exist for Equation (H.24) are 


*3 


» 0 


(H.21) 


(H.22) 


(H.23) 


(H.24) 


(H.25) 
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*03 " *02 * ® 


(H.26) 


<j ) 3 - »* 


(H.27) 


where n = ± 0,1, 2,3 Avoiding the trivial solutions, pick n - - 2 so that 


4>, - -2it. 


(H.28) 


From Equation (H.8), the third impulse occurs in time at 


2 n 


‘03 




(H.29) 


Substituting Equations (H.22) and (H.28) into Equation (H.10) yields 


fi 2 sm(4> 2 ) * 0. 


(H.30) 


The two solutions that exist for Equation (H.30) are 


B 2 - 0 


(H.31) 
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<j> 2 ■ Wit 


(H.32) 


where n = ± 0, 1,2,3,... . Avoiding the trivial solution, pick n - -1 so that 

4» 2 ■ - « . (H.33) 


From Equation (H.8), the second impulse occurs in time at 


u.l/W 5 


(H.34) 


Now that the times at which the impulses occur is known, the amplitudes of the four 
impulses can be found. Substituting Equations (H.22), (H.28) and (H.33) into Equations 
(H.9), (H. 11) and (H.13) produces 


B l -B 2 +B 3 -B 4 - 0 


(H.35) 


~^2 r 02 *B 3 t Q3 -B 4 t 04 * 0 


(H.36) 


” 82*02 + 83*02 ~ 84*04 * 0 • 


(H.37) 
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■ 


From Equations (H.23), (H.29) and (H.34), the relationships between the times are 


*0 3 “ 2*02 


(H.38) 


t 


04 


3r 


02 ' 


(H.39) 


Substituting these time relationships into Equations (H.36) and (H.37) yields 


B 2 « 3Bj 


(H.40) 


B 3 - 3 B x 


(H.41) 


*4 



(H.42) 

* 


At t = Equation (H.40) produces 


JTT 2 v'w 5 


(H.43) 
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or 


c» 

A 2 - 3 A x e '^ rr? . 


(H.44) 


At l - fo,, Equation (H.41) produces 


/w 3 v'w 5 


(H.45) 


or 


ac« 

A 3 - 3A,e 


(H.46) 


At / = t^. Equation (H.42) produces 


m - C«. - » M ) = _ - C«. ('m - »•!> 

/w 5 /I 7 ? 


(H.47) 


or 


^4 


A,e 


3C« 

V'l-C* 


(H.48) 
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If the amplitudes are normalized so that 


£4- 1 (H- 49 ) 


and let 


M 



(H.50) 


then the amplitudes of the four impulses and the times at which they occur are 




i (H.51) 

l+3M + 3M 2 +M 3 


t n « 0 (H.52) 


*2 


3M 

1 +3M + 3Af 2 +Af 3 


(H.53) 


*02 


It 

o./W 5 


(H.54) 


^3 - 


3M 2 

1 +3M + 3Af 2 +Af 3 


(H.55) 


2n 


*03 


-c J 


(H.56) 


^ = 


Af 3 


4 1 + 3Af + 3Af 2 + M 3 


(H.57) 


3 it 


*04 




(H.58) 
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Now that the amplitudes and times of the impulses are known, the theoretical impulse 
response of a second-order system can be found using the results from Appendix D. 
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APPENDIX I 


CALCULATIONS FOR DEVIATIONS IN NATURAL FREQUENCY 

This appendix evaluates the residual vibration of the input shaping technique when 
the actual natural frequency, u a , deviates from the design natural frequency, The 
magnitude of the residual vibration is important since the impulse characteristics, used 
for input shaping, are based on the design natural frequency values. 

Using the results derived in Appendix D, the actual response of a second-order 
system to an impulse of unity gain is 

x a (t) = — - C 2 0- fl- 1 ) 


If the single impulse is divided into two impulses with characteristics derived in 
Appendix F, the actual response for t ^ t m becomes 


x 2a (0 - x l (t)*x 2 (t) 


( 1 . 2 ) 
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where 


*,(f) 


1 

Af + 1 




/w 5 


«' Cw ‘*sm(o» a v^l - C 2 f ) 


a-3) 


x 2 (r) 


w - 

Af + 1 * 

/R 


C w # (i - ipj) 




(1.4) 


Using the results derived in Appendix E, the actual system response to the two impulse 
input for * ^ *02 can be expressed as 


* 2 ,(0 “ B 2a sin(at+y 24t ) 


a.5> 


where 


B 


2 a 


1 

(i) 

Af + 1 


/w 5 


M 


1 +2Af 





a-6) 
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To determine the amount of residual vibration, a vibration error expression must 
be defined. The error, denoted err , can be expressed as the ratio of the actual multiple 
impulse response magnitude to the actual impulse response magnitude of a second-order 
system. The error expression is defined only for time after the multiple impulse input 
has occurred to ensure that the system has received identical amplitude inputs. 
Mathematically, this is written as 

M . for t ± , ot M> 

|*.(<)| 

where k is the number of impulses. The residual vibration is the vibration error 
expressed as a percentage. Using Equations (1. 1) and (1.5), the vibration error for the 
two impulse input becomes 
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If the single impulse is divided into three impulses with characteristics derived in 


Appendix G, the actual response for * ^ *03 becomes 

*,.(<) * X,(t) ♦*!(») .*j(») 


a. io) 


where 


x,(r) 


1 


1 *2 M + M 2 


G) 


/w 5 


Cw * r sin (« # Vl ” C 2 t) 


(I.H) 


x.(r) 


2M 

1 + 2M + M * 

/w 3 


u. 


C«.(r- 


<f 1 - C 2 (f - f 0 i)) 


(1.12) 


A/ 2 


xj(o - a.i3) 

/w 5 


Using the results derived in Appendix E, the actual system response to the three impulse 
input for t 2: can be expressed as 


* 3 ,(0 * B ia sm(at+^f u ) 


a. i4) 
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where 


(»> 


B 


3a 


1 + 2Af + M 1 m -C 

’ "" v 




a. is) 



a- 16) 


with 


nmn 3- 



a.i7) 


^»3a " 


1 + 2Af 


K) 






(U8) 


Substituting Equations (1.1) and (1.14) into Equation (1.8), the vibration error for the 
three impulse input becomes 


|*.(0| 


1 

1 + 2Af + Af 2 


>/ ( ) 2 + ( ) 2 ‘ 


(1.19) 
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If the single impulse is divided into four impulses with characteristics derived in 


Appendix H, the actual response for r ^ r M becomes 

*4,(0 “ * i (0 + * 2 (0 + * 3 (*> +X 4<*> 


where 


*,(*> 


1 

1 + 3M + 3Af 2 + M 3 

/W 5 




sm^o^/T-C 5 r ) 


x 2 (r) 


3M 

V i-c 1 






-C 2 (f - f 02 )) 


3A/ 2 


<0. 


*,(0 - 1 

/w 5 






x 4 (r) = 1 1 3Af 1 3 M 2 + M 3 e ~ c ^ (, ~ <t 4 ) siii(co a Vl - C 2 (* - r 04 )) 

/W 5 


a.20) 


a.2i) 


( 1 . 22 ) 


(1.23) 


(1.24) 
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Using the results derived in Appendix E, the actual system response to the four impulse 
input for r ^ can be expressed as 


* 4.(0 “ 


a. 25) 


where 


B 


Ac 


l 

1 + 3Af + 3Af 2 + Af 3 

/T^c 5 




' C "‘ 7 ( An 4-) 2+ ("“ W <-) 2 


a. 26 ) 


- “tea’ 1 


' \ 

V *”4. j 


a.27) 


with 






(1.28) 
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Substituting Equations (1.1) and (1.25) into Equation (1.8), the vibration error for 


the four impulse input becomes 


* 4«(01 

|*.CO| 


1 

1 +3A/ + 3A/ 2 


\/( <feB 4 -) 2 + («“ W 4.) 2 * 


(1.30) 
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APPENDIX J 


CALCULATIONS FOR DEVIATIONS IN DAMPING RATIO 

This appendix evaluates the residual vibration of the input shaping technique when 
the actual damping ratio, deviates from the design damping ratio, f. The magnitude 
of the residual vibration is important since the impulse characteristics, used for input 
shaping, are based on the design damping ratio values. 

Using the results derived in Appendix D, the actual response of a second-order 
system to an impulse of unity gain is 

x„(0 - “* « " t ) . (J.l) 

If the single impulse is divided into two impulses with characteristics derived in 
Appendix F, the actual response for t ^ t m becomes 

x 2fl (f) * x i( f ) +x i(0 (J 2) 
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where 




(J.4) 


Using the results derived in Appendix E, the actual system response to the two impulse 
input for t ^ *02 can be expressed as 

“ fi 2a sin(«f+t^) (J.5) 

where 
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To determine the amount of residual vibration, a vibration error expression must 
be defined. The error, denoted err , can be expressed as the ratio of the actual multiple 
impulse response magnitude to the actual impulse response magnitude of a second-order 
system. The error expression is defined only for time after the multiple impulse input 
has occurred to ensure that the system has received identical amplitude inputs. 
Mathematically, this is written as 


err 



for f a t ok 


(J.8) 


where k is the number of impulses. The residual vibration is the vibration error 
expressed as a percentage. Using Equations (J.l) and (J.5), the vibration error for the 
two impulse input becomes 



(1-9) 
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If the single impulse is divided into three impulses with characteristics derived in 


Appendix G, the actual response for becomes 

x 34l (0 - Xj(f) +x a (0 +x 3 (f) 


(J.10) 


where 


x,(f) 


1 

1 ♦ 2Af + Af 2 





-fjf) 


a. id 


2Af 


0 ). 
r2 * 


/w! 


x,(0 « 


Af 2 

1 + 2Af + Af 2 




-£(f - f 03 )) . (L13) 


Using the results derived in Appendix E, the actual system response to the three impulse 
input for r ^ f(u can be expressed as 


x 3a (0 - B 3a sm(af+* 3 ,) 


d. 14) 
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where 


B 


3 a 


— e - W J(den u f + (num 3a ) 2 

yi-c 


a. i5) 



rnw^A 

3« ) 


(J.16) 


with 


rtum 3a * 2Af 


(i.i) | 

\ « /sin 


l*N 


i-ci 

T T>) 


+ M 


HL 


2 ic 


1-C. 
N i-c 2 J 


a. i7) 


den 


3 a 


(-T> 

= 1 +2M' c 'cos 



^ 

i-ci 

2 ( 

•4) 

r 

\ 

i-cS 

It 

N 

1-C 2 J 

+ Af ' 

€ 7 cos 

2 it 

N 

i-c 2 j 


(J. 18) 


Substituting Equations (J.l) and (J.14) into Equation (1.8), the vibration error for the 
three impulse input becomes 


|*.<0| 


1 

1 +2M+M 2 


v f(den 3- ) 1 +(nu , n 3- ) 2 . 


(J.19) 
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If the single impulse is divided into four impulses with characteristics derived in 


Appendix H, the actual response for t becomes 

*4«(0 - VO + VO + VO + VO ( J - 2 °) 


where 


«,(0 . 0 . 2 D 

ft 7 ?. 


VO 


3M 

1 +3M + SM 2 +M 3 


<l> 
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Using the results derived in Appendix E, the actual system response to the four impulse 
input for t ^ can be expressed as 


where 


x 4 ,(0 “ B 4a sm(af+t 4- ) 
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Substituting Equations (J. 1) and (J.25) into Equation (J.8), the vibration error for 
the four impulse input becomes 


|*,( 0 | 
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1 + 3Af + 3M 2 +M 3 


/(*»4,) 2+ ( n “ m 4,) 2 • 


(J.30) 
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APPENDIX K 


CIRCUIT FOR ACQUIRING FREQUENCY DATA 


The circuit shown in Figure K. 1 was used to combine the control signal from the 
MicroVAX with a random noise signal from the HP Analyzer. This new signal was then 
input to the KEPCO amplifier that drives the hydraulic valve for joint one. 
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Figure K. 1 Circuit Used to Combine Random Noise with Control Signal 


The resistor values are chosen to insure proper impedance matching of the various 

system components. The resistor values used are 

R, = 10 kG 
R 2 = 0.1 - 10 kO. 

NOTE: The KEPCO amplifier requires its input device to have an impedance higher than 
the operational amplifier (op amp) can supply. Therefore, R 2 must be used to increase 
the overall back impedance of the op amp. 
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APPENDIX L 


LEAST-SQUARES REGRESSION OF EXPERIMENTAL DATA 


This appendix shows how a least-squares curve fit was derived from the 
experimental data. A polynomial regression that is a function of B x and d 2 will be 
developed for the damped natural frequency, f d , and the damping ratio, f. A full 
discussion of linear regressions can be found in [10]. 

Consider the following S^-order polynomial 

y * bjOj + + + + fl o • (LI) 


The sum of the squares of the residuals is 


~ a 0 ~^1®2 i ~ ^2®2 i ~ ^3 ®2i) 


(L.2) 


where N is the number of data points. Take the derivative of Equation (L.2) with respect 
to each of the unknown coefficients of the polynomial to get 
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H 


Equations (L.3) through (L.9) can be set equal to zero to minimize S r Rearranging the 
previously mentioned equations yields a linear equation with respect to the unknown 
coefficients of the form 

[A]Z - z ( L * 1 °) 

where 



N 

£»„ 

i-l 

N 

E e ii 

i-l 

AT 

E 0 ?. 

i-l 

N 

£9 31 

i-l 

Af 

E 0 i 

i-l 

CD 


N 

£9 „ 

i-l 

N 

E el 

i-l 

E«t 

i-l 

Af 

E 0 fi 

i-l 

£91,92. 

i-l 

ffl.,91, 

i-l 

£ 01.92. 

i-l 


S 

E ef, 

i-l 

Eel 

i-l 

E e f< 

i-l 

AT 

E 0 u 

i-l 

£e»92, 

i-l 

E®li 0 2i 

i-l 

E ®?|02i 

i-l 

[A] = 

M 

E e u 

i-l 

E«u 

i-l 

AT 

E 0 i 5 < 

i-l 

N 

E 0 t< 

i-l 

£e?,e 2 , 

i-l 

E 0 li 0 2< 

i-l 

E 0 W 0 2< 

i-l 


N 

E ®2i 

i-l 

E 0 2i 0 li 

i-l 

i-l 

£9*9?, 

i-l 

N 

E 0 2 2 i 

i-l 

N 

E 0 2i 

i-l 

E 0 i 

i-l 


N 

Eel 

i-l 

E 0 2< 0 li 

i-l 

Zelef, 

i-l 

E 0 2i 0 ?i 

i-l 

N 

E 0 2i 

i-l 

N 

E 0 2 4 i 

i-l 

N 

E 0 2< 

i-l 


N 

E®1 

i-l 

E ®2i®li 

i-l 

E 0 2i 0 fl 

i-l 

E e 2 3 i e ii 

i-l 

N 

E 0 i 

i-l 

H 

E 0 i 

i-l 

E 0 2* 

i-l 


CL. 11) 


139 




140 


The unknown coefficients can now be easily solved for by 


Z-[AV l z. < L - 14 ) 

A MATLAB program was used to evaluate Equation (L.14) using the 
experimental data from Chapter IV. The resulting equations for the damped natural 
frequency and damping ratio are 

//0 lt 0 2 ) = - 0.4212 0 2 + 1.8086 6j - 2.1441 0 2 - 0.222901 ^ J5) 

+ 1.4540 0j - 3.6392 0j + 7.7830 


and 


C(O lf 0 2 ) = - 0.4269 0 2 + 1.8918 0 2 -2.71590 2 -0.053301 
+ 0.3410 0j - 0.6148©! + 1.6762 . 


(L.16) 
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